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INTRODUCTION 


Development  of  nitramines  as  solid  propellant  ingredients  has  identified 
the  sketchy  understanding  of  their  combustion  mechanisms.  The  chemistry  and 
hydrodynamics  have  yet  to  be  satisfactorily  combined  to  give  a  description  of 
how  the  combustion  process  could  be  controlled  or  modified.  Present  develop¬ 
ments  rely  on  trial-and-error  empiricism  in  adopting  formulations  and  predict¬ 
ing  their  response  to  variations  in  pressure,  initial  temperature,  gas  velocity, 
and  composition.  Improvements  are  needed  in  the  theoretical  description  of  the 
combustion  process.  A  model,  developed  at  Princeton  University,  forms  the 
basis  for  improvements. 

The  objective,  in  general,  is  to  improve  the  theoretical  description  of 
nitramine  propellant  combustion,  through  development  of  a  comprehensive  model. 
More  precisely,  this  is  to  be  accomplished  through  the  following  tasks: 

(1)  Devise  an  improved  physical  model  for  the  vaporization  and  liquid 
layer  of  a  nitramine  monopropellant.  Use  the  improved  model  to  compute 
regression  rate  dependence  on  pressure  and  initial  temperature  over  the 
pressure  range  7-70  MPa  and  temperature  range  230-330  K. 

(2)  Devise  an  improved  submodel  for  the  liquid  layer  accounting  for 
bubbles.  Use  the  improved  model  to  calculate  combustion  solutions  to  be 
compared  with  earlier  published  solutions  at  1-^4  MPa. 

(3)  Perform  a  critical  review  of  other  combustion  models  of  monopropellant 
and  of  nitramine  based  composite  propellant. 

(k)  Perform  a  critical  analysis  of  the  effect  of  the  three  chemical 
reactions  used  in  the  existing  monopropellant  model  with  a  view  to  a 
scheme  to  catalyze  the  reactions. 

(5)  Collect  and  analyze  available  data  on  nitramine  monopropellant 
combustion  with  a  view  to  suggesting  which  data  should  be  considered 
baseline  for  model  validation. 

Items  (3)  and  (5)  are  performed  in  Section  2,  and  items  (l)  and  (2)  in 
Section  (3)  herein.  Task  (U)  is  still  incomplete  at  the  time  of  this  writing;  an 
effort  is  made  to  conclude  this  study  in  the  near  future. 


I.  CRITIQUE  OF  THE  BECKSTEAD-DERR-PRICE  (BDP) 
PROPELLANT  COMBUSTION  MODEL 


1 . 1  Background 

Over  a  decade  ago,  Becksteaa,  Derr  and  C.F.  Price  introduced 
a  model  of  solid  monopropellant  and  propellant  combustion  £  1 , 2] , 
referred  to  as  BDP  herein.  The  corresponding  flame  structure 
postulated  is  depicted  in  Fig.  1,  where  the  overall  physical  concept 
is  shown.  It  consists,  for  each  oxidizer  crystal  exposed  at  the 
surface,  of  three  types  of  gaseous  flames.  As  will  be  discussed 
later,  each  of  these  flames  is  assumed  to  contribute  a  component  of 
heat  feedback  proportional  to  exp  (-z .  )  where  z.  (i  =  1,  2,  3) 
denote  dimensionless  Mflame  standoff"  distances.  The  BDP  approach 
was  used  widely  in  the  U.S.  ever  since,  and  was  applied  to  a  rather 
diverse  propellant  lot,  including  ammonium  perchlorate  composites, 
composite-modified  double  base,  and  nitramine  {KMX ,  RDX)  propellants. 
Interestingly,  a  previous  composite  propellant  combustion  model,  the 
granular  diffusion  flame  (GDF)  model  of  Summerfield  [3,4,53  was  far 
less  accepted  in  comparison,  despite  its  considerable  success  in 
simulation  of  burning  rate  pressure-dependence  for  a  large  class  of 
ammonium  perchlorate  composites  163.  The  GDF  model  was  criticized 
for  its  intuitive  formulation,  [7,8j  based  on  dimensional  or  physical 
scaling  arguments,  rather  than  formal  derivation;  its  success  in 
burning  rate  simulation  was  attributed  to  the  availability  of  two 
adjustable  constants. 

An  appreciable  number  of  well-developed  composite  propellant 
combustion  models  has  been  reported  over  the  past  decade,  which 
endeavored  to  incorporate  many  real-life  effects,  such  as  propellant 
structure,  energetic  binder  and  aluminized  configurations.  Such  are 
the  works  of  Beckstead  [9-13],  Cohen  [14-17],  King  [18,19],  and  the 
statistical  refinements  by  Glick  [  20-23]  ,  all  of  which  had  evolved 
from  the  original  BDP  model  of  1970.  In  most  cases,  good  comparisons 
with  experimental  burning  rate  vs  pressure  were  obtained. 
Surprisingly,  the  gas  phase  combustion  effects  are  still  in  the 
exp  (-z^)  form. 


In  contrast  to  the  highly  approximate  treatment  of  combustion 
processes  in  all  of  the  BDP-type  models,  the  following  aspects  assume 
great  relative  importance:  propellant  matrix  structure,  oxidizer 
particle  sizes  and  distributions,  and  surface  structure  (e.g.,  convex 
or  concave  burning  oxidizer  surfaces,  appearance  of  melt  layers) . 

This  is  not  to  say,  unfortunately,  that  these  structural  eleirients 
have  been  rigorously  treated.  To  the  contrary,  their  representation 
and  relevance  to  the  burning  process  can  be  inferred  only 
speculatively  at  present.  Indeed,  the  incorporation  of  these 
properties  gave  rise  to  a  large  number  of  adjustable  constants  in 
the  formulation,  which  can  not  be  derived  analytically  nor  verified 
by  physical  measurements. 
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To  obtain  some  preliminary  insight  into  the  development  of  BDP 
modeling  over  the  years,  the  following  calculation  is  carried  out, 
cum  granu  salis.  Regarding  the  surface  energy  balance  (the 
centerpiece  of  all  BDP  models) ,  the  symbols  denoting  combustion 
effects  (i.e.,  the  exponential  terms),  and  the  totality  of  right-hand 
side  symbols  (gas  phase  heat  feedback)  are  enumerated;  their  ratio 
(the  so-called  combustion  fraction)  represents  a  point  in  Fig.  2  for 
each  investigation.  The  resultant  decaying  linear  trend  obtained 
by  least-squares  fit  demonstrates  that  the  relative  importance  of 
combustion  in  solid  propellant  regression  is  diminishing  in  time  and 
may  even  vanish  during  the  1980s.  One  may  envision  the  Space  Shuttle 
boost  phase  in  1990  powered  entirely  by  oxidizer  granularity,  modality 
and  volume  fraction. 

The  foregoing  trend  should  be  envied  by  other  combustion  system 
models,  but  deserves  closer  examination:  Is  the  exponential 
combustion  term,  as  stated  by  BDP,  consistent  with  modern  combustion 
theory,  available  since  the  1950s?  Are  the  physical  and  analytical 
BDP  models  compatible  with  their  own  physical  picture,  as  depicted 
in  Fig.  3?  These  are  some  of  the  questions  which  motivate  the  present 
critique.  More  precisely,  regarding  Fig.  3,  the  questions  posed  for 
this  analysis  are:  (1)  Is  the  monopropellant  flare  properly 
formulated  mathematically?  (2)  Can  the  heat  feedback  contributions 
from  the  various  "adiabatic"  flame  components  be  simply  superimposed, 
without  interaction?  (3)  Can  a  planar  flame  sheet  at  a  mean  standoff 
distance  represent  a  conical  or  parabolic  diffusion  flame  surface? 

(4)  Is  this  parabolic  flame  surface  truly  diffusion  controlled? 

The  approach  is  simply  to  check  the  validity  of  the  BDP  combustion 
model  within  the  framework  of  its  own  physical  picture  or  basic 
assumptions.  This  is  carried  out  in  a  constructive  manner,  with 
parallel  derivations  offered  which  are  in  agreement  with  present 
knowledge  of  combustion  theory  and  implement  available  thermochemical 
data. 
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1 , 2  The  Monopropellant  (Oxidizer)  Flame 


The  basis  for  the  heat  feedback  expression  in  the  BDP 
irtonopropellant  analysis  C  1 , 2] ,  viz, 

d'C /cte(o+)  -  €T2b  (1) 

is  an  approximation  of  the  single  (overall)  exothermic  chemical 
reaction  step  in  the  gas  phase  by  a  Dirac  delta  function,  at  a 
dimensionless  standoff  distance 


%=  ^/(X/wCp)  (2) 


The  energy  interface  condition  (at  z  =  0+)  corresponding  to  Eq.  (1) 
as  used  by  BDP  is 


dr/cte (o+)  =  [cc (rs -To )~Qt  ]/Q*?  (3) 

where  the  dimensionless  thermal  enthalpy  is  C  =  C  (T-Tc)/Q*,  C*  is 
the  net  surface  heat  release,  and  Q*  is  the  net  g§s  phase 
exothermicity .  The  result  of  Eq.  (1)  obviously  pertains  to  the 
solution  of  the  homogeneous  energy  equation  =  0  for  the  region 

0  <  z  <  z^ .  This  strictly  precludes  the  effect  of  chemical  reaction 
and  necessitates  extraneous  introduction  of  chemical  kinetics  effects 
through  the  intuitive  standoff  approximation  y*  *  m /XI  f  proposed 
earlier  by  Hermance  C24].  Thus,  Eq.  (2)  yields 

zB  =  (4> 

where  m  is  the  mass  flux  and  JT2  denotes  the  mean  reaction  rate. 


Alternative  derivation  (KPR) 


The  use  of  this  extraneous  condition  can  be  entirely  avoided 
when  a  more  rigorous  analysis  (still  within  the  delta  function 
approximation  framework)  is  carried  out,  as  follows.  Consider  the 
following  simplified  chemical  reaction  scheme: 


SUBSURFACE  DECOMPOSITION: 

R(c) 

— *■ 

P(g) 

+ 

Qc 

(5. a) 

VAPORIZATION  @  SURFACE: 

R(c) 

R(g) 

+ 

(5 .  b) 

GAS  PHASE  REACTION 

R(P) 

P(g) 

+ 

(5.c) 
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R  and  P  denote  reactant  and  product  mixtures, respectively. 

0  *  Q,  +  Q  is  the  net  heat  of  reaction  in  the  condensed  phase, 

c  1  v 

and  vaporization  in  step  (5b)  can  be  either  simple  phase  change  or 
dissociative.  Note:  Q>0  denotes  exothermicity  herein.  The  physical 
model  for  the  delta  function  flame  is  shown  in  Fig.  4. 

The  dimensionless  conservation  equations  in  the  gas  phase,  for 
0  <  z  <  co  are 


(6-a> 

dV/bte  =  -ASO-^)  (6*b) 

where  0  =  C  (T-T  ) /Q*  and  Y  denote  reduced  thermal  enthalpy  and  a 
reference  (consumable)  species  mass  fraction,  respectively.  The  flame 
speed  eigenvalue  is  defined: 

A  5  ^  JQGfjTp)/™2-  (?) 


with  -H  being  the  mean  reaction  rate.  The  boundary  data  are 


e  co+) = o ,  de/<te«>+)  =  <j,0 ,  6(«) = e4  <e ' 

YC0+-)  -dv/d2fo+)  =  6  (9-a| 

't'(oo)  =  dV/c(2.(o°)  =  0  .  <9.M 

where  0^  =  Cp(Tf  -TsYQ*  ^ 

%=l  CcCrs -t0 )  -OM /q; 


Using  an  integrating  factor,  the  solution  of  the  differential  system, 
Eqs.  (f)-(9)  is  readily  found 


0(e)  =  1 )  -  A(e?  1)  <u-a> 

where  1  (z-z^)  is  the  Heaviside  (or  unit  step)  function?  a  similar 
solution  obtains  for  the  species  equation.  A  single  integration  of 
Eq.  (6 . b)  between  0  and  infinity  is  facilitated  by  the  delta  function 
imposed.  This  yields,  after  implementation  of  Eqs.  (9.a,b): 

A.  =  (11. b) 
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where  G  is  the  total  incoming  flux  of  reactant  'R,'  available  at  the 
surface.  Note  that  G  <1  whenever  subsurface  decomposition  prevails. 
The  net  surface  heat  release  is,  cf  Eqs.  (5a,b,c): 


Q*  ~  —  +-  (12) 


Note  that  all  subsurface  heat  release  assumed  to  occur  at  the  surface, 
and  the  total  gas  phase  heat  release  is 


(13) 


In  general,  the  extent  of  subsurface  decomposition  depends  on  surface 
temperature,  residence  time  and  other  parameters:  G  =  G  (T  ,  m;  T  ) 
can  be  written.  Hence  G  may  be  variable;  this  has  not  been£recognized 
in  the  EDP  model,  which  precludes  any  consideration  of  chemical 
species. 

The  heat  feedback,  obtained  from  a  single  integral  of  Eq.  (6. a), 
using  an  integrating  factor: 

cis/d^  (14) 

is  formally  similar  to  that  of  BDP ,  cf  Eq.  (1).  However,  the  central 
issue  of  the  present  derivation  as  regards  burning  rate  calculation, 
is  the  flame  speed  eigenvalue, 

jflC'PjTf  )//W2-  =  G>  (is) 


combining  Eqs.  (7)  and  (11),  which  naturally  incorporates  the 
subsurface  processes  along  with  the  effect  cf  gas  phase  chemical 
kinetics. 

The  mean  reaction  rate  in  Ea.  (15)  can  be  put  more  explicitly, 

-QCp/Tl)  =  as) 

where  and  (3^  denote  the  overall  reaction  apparent  order  and 
activation  energy,  respectively,  while  Kt  is  proportional  to  the 
prefactor.  Substitution  in  Eq.  (15)  yields 
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The  BDP  analysis  obtains  an  analogous  expression,  using  the  heat 
feedback  expression,  Eqs.  (1),  (3)  and  (4): 


Comparison  between  the  burning  rate  expressions,  Eqs.  (17,18)* 
reveals  a  fundamental  difference,  demonstrated  by  the  following  simple 
test  case.  Whenever  G (T  ,  m;  T^)  =  const,  and  =  const,  the  region 
0  <z  <oo  is  equivalent  lo  the  same  region  in  a  pure  gaseous 
deflagration  (subject  to  the  same  delta  function  approximation,  of 
course) ,  regardless  of  "surface  temperature"  variation. 

In  this  instance,  m^pn^2  is  expected,  as  indeed  qiven  by  Eq.  (17) 
derived  herein,  but  not  from  the  analogous  BDP  expression,  Eq.  (18), 
where  an  explicit  T  -dependence  prevails;  this  dependence  indicates 
inconsistency  with  the  physics  of  the  problem. 

Note  that  the  validity  of  the  foregoing  equivalence  argument 
(with  pure  qaseous  def laqration)  can  be  readily  established 
mathematically  by  noting  that  G  =  const  implies  (letting 
C  =  C  for  this  purpose) 


c  p 


(19) 


as  obtained  by  integration  of  Eq.  (6. a)  once,  and  using  Eqs.  (8);  this 
is  compatible  with  the  overall  enthalpy  balance.  Consequently,  a 
manifold  of  genuine  solutions  can  still  be  generated  for  distinct 
values  of  ,  as  long  as  the  linear  combination  of  Eq.  (19)  is 
maintained,  each  solution  with  a  distinct  burning  rate  (or  flame 
speed) ,  despite  the  eigenvalue  A-  being  fixed. 

It  should  be  pointed  out  that  the  delta  function  flame  treatment 
herein  is  quite  crude,  and  was  intended  only  to  simulate  the  physical 
model  used  by  BDP,  Contrary  to  popular  legend,  propagated  by  an 
erroneous  footnote  in  the  KTSS  paper  £253,  this  approach  was  never 
initiated  by  von  Karman;*  it  was  used  by  Zeldovich  in  a  simplified 
treatment  of  chain  reactions  [27].  There  exist  much  more  rigorous 
analytical  treatments  of  gaseous  def lagrations ,  using  matched 
asymptotic  expansions,  by  Jain  and  Kumar  £28],  Bush  and  Fendell  [29], 
Joulin  and  Clavin  [30]  and  Ben-Reuven  [31]. 

To  obtain  formal  closure,  by  which  the  burning  rate,  m,  could 
be  calculated,  both  formulations  require  a  pyrolysis  law,  or 
m  =  m  (Tc) ;  in  addition,  Eq.  (17)  requires  also  explicit  knowledge 

of  subsurface  processes,  to  obtain  explicitly  G  (T  ,  m;  T  ) .  Once 

s  o 


*  in  particular,  no  mention  of  such  approximation  was  made  in 
the  classical  von  Karman  paper  given  at  the  Sixth  International 
Combustion  Symposium  in  1956  [26]. 
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closure  is  obtained,  the  burning  rate  pressure  sensitivity  and  the 
temperature  sensitivity  can  be  derived  for  both  models.  In  general, 
these  properties  are  defined. 


W/™  I 

11  =  Sr/  Mt0  j 


j 

STo 


(20) 


where  ^  here  denotes  the  variational  operator. 


From  Eq.  (17),  the  first  variation  leads  to 


S&  ,  7 
&  w 


(21) 


where 


Cp<' 


-£'/Tf 


K,  =  A1T“,'/(RuT)’ni  (22) 


Note  that  we  assume  that  the  product  (  A/C  )K4  is  independent 
of  temperature  for  any  reaction  order  n^  un8er  consideration.  Now, 

Scb  __  r  21L  +  ”|  Sm  <-  (23  . 

*3"  “  L  3Ts  J 

Sk:  __  r  ^  ^Ts  "1  ilf  \r  (23. b) 

If  "  L  9Tf  dTs  BiuvnJ  VK  3Tf  dX0l° 

The  algebraic  sum  of  all  Sm/m- terms ,  collected: 

—  dTs  ,  'dL>J=^  _  5ZL.92L  +9  (23. c) 

^  “  'a-T^  ■*“  *  c>7* 

Thus,  using  Eqs.  (20),  (21)  and  (23): 

(24*a) 

cfpCMBR)  =  (  H  “  W^A  (2,-b| 

which  still  require  explicit  definition  of  G  (T  ,  m?  T  ) .  In 
comparison,  the  corresponding  BDP-derived  sensitivities  can  be  written 
explicitly, 


n(BDP)  =  U,  !\z  +  crTsY(ilQ{  e6€  J 
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$CBaP)  =  MtlSMzdl  . 

Y  z  +  q»T <25-b> 

where,  cf  Eqs.  (1) , (3) , 

^  r  CcC-ll-To)-C)r  j 

Note  that  C  =  C  is  assumed  in  the  original  BDP  formulation,  whereby 
is  formally  independent  of  T  . 

Evidently,  the  actual  values  obtained  by  Eqs.  (24. a, b)  herein 
depend  on  the  particular  subsurface  processes,  whereby 
G  =  G  (T  ,  m;  T  )  is  defined,  which  is  expected  to  differ  from  one 
monopropellant  type  to  another.  In  contrast,  the  sensitivities 
obtained  from  the  BDP  formulation,  Eqs.  (25. a, b),  depend  only  on  a 
suitable  choice  of  the  parameter  family,  (E^ ,  Q* ,  Q*,...),  while  the 
underlying  functional  dependences  remain  fixed  for  all  monopropellant 
types.  This  is  a  rather  strong  statement,  regarding  the  universality 
of  monopropellant  decomposition  and  combustion  mechanisms,  totally 
unwarranted  by  the  physical  processes  involved,  '  The  foregoing 
statement  does  not  detract,  however,  from  the  simulative  capabilities 
of  the  function 

-PB(XS)  =  Ti  v eg  e  (26) 

which  appears  in  the  denominator  of  Eqs.  (25. a, b).  As  will  be 

demonstrated  later,  f  has  a  remarkable  correlative  power. 

B 

To  facilitate  comparison  with  the  BDP  model,  two  highly 
simplifed  cases  are  constructed,  as  follows. 

(1)  For  AP  simulation, 

GCTs'M  )=  0, 3  =  Const,  7^  =  500  <  (27. a) 

i.e.,  70%  of  AP  decomposition  is  assumed  to  occur  with  the  subsurface 
region,  in  line  with  similar  considerations  by  BDP,  and,  in  a  much 
more  detailed  AP  deflagration  model,  by  Guirao  and  v;illiams  £32]. 

In  addition,  C  >  C  is  assumed  herein,  leading  to  an  explicit 
dependence  of  ?he  aHiabatic  flame  temperature,  T^,  on  T^  through 
the  overall  enthalpy  balance: 

CcOTs~lo)  +  £f(Tf  =  G^Qsj  +  (27'b) 

whereby, 

=  (Cp  -Cc)/  cp  (27.0 

Note  that  although  (C^  -  C  )T  /q  is  expected  to  be  small,  still 

(C  -  C  )/C  ~  0(1)  accord?nqSto  available  data, 
c  p  p 
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It  should  be  pointed  out  that  $G  =  0  at  =  0,  while  ^Tg, 

£m,  £>p  £  0,  may  still  follow  from  the  general  case  outlined  in 
Eg.  (23 .a) : 

'dJU<£> .  _  n 

'dS^yr\ 

This  implies  that  residence  time  and  characteristic  chemical  kinetics 
relaxation  time  in  the  subsurface  region  balance  exactly  to  yield 
zero  variation  of  G  under  these  conditions.  However,  for  £Tq  f  0, 
while  =  0  and  £T£,  £m  ^  0,  one  may  have 

■SG  = 

<3  0T& 

from  Eq.  (23. a).  The  corresponding  correction  to  the  temperature 
sensitivity  in  this  particular  case  is  derived  in  the  second  part 
of  Appendix  A. 

(2)  For  nitramine  simulation  (HMX  in  particular) ,  the  extent  of 
subsurface  decomposition  expected  is  much  smaller  in  comparison,  but 
variable.  The  following  approximation  is  based  on  the  assumption 
that  the  activation  energy,  F  ,  is  high,  and  hence  most  of  the  first- 
order  overall  reaction  in  the  liquid  phase  occurs  in  a  thin  region 
near  the  surface.  For  the  derivation  given  in  Appendix  A,  the  result 
is 

<SC^,miTt,)±e>Cp[-kc'€l/e£ 

Kc=  Ac?cAt  /Ccmz  ,  o*£=Ts/pc«  i 

The  appropriate  partial  derivatives  are 

[  £  +  IjpJ  ]  A 

'oh^(o/'drf\  —  —'2.1/^Gj /yr\ 

rc  =  A6/(ts-t„) 

Further,  in  this  case,  C  ==  C  is  considered  in  good  agreement  with 

measured  C  and  calculated  C  ^(gas  phase  near-field)  data, 
c  p 

In  both  cases  (1)  and  (2),  the  following  pyrolysis  law  is 
imposed,  to  obtain  closure: 


(28. a) 

(28. b) 

(28. c) 


ynfji)  =Ase 


(29) 
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One  may  proceed  now  to  compare  the  presen*-  derivation  and  the  BDP 
model,  for  AP  and  HMX;  the  appropriate  input  data  are  given  in  Tables 
1  and  2.  Experimental  burning  rate  data  for  HMX  and  AP  are  taken 
from  Price,  Boggs  and  Derr  t 333,  and  for  AP  also  from  Hightower  [343. 


TABLE  1 


AP  INPUT  DATA 


PROPERTY 

UNITS# 

PRESENT 

MODEL 

BDP  MODEL, 

PRESENT  CALCULATION 

\ 

C 

cal/cm-s-K 

9xl0“4  (c) 

- 

C 

C 

cal/g-K 

0.5  (b) 

0.3 

Pc 

g/cc 

1.94  (a) 

- 

A 

c 

1/s 

2.39xl018  (c) 

- 

E 

c 

kcal/mol 

58.22  (c) 

- 

n 

_ 

0  (c) 

_ 

c 

As 

kg/m^-s 

1 . 29xl07  (d) 

3xl06 

E 

s 

kcal/mol 

24.24  (d) 

22 

A  /C 

9  p 

kg/m-s 

5x10  3  (e) 

10~3 

C 

P 

cal/g-K 

0.3  (c) 

0 . 3=C 

c 

°c=vQiub 

cal/ g 

+444.6  (d) 

Q*=+210 

Q!ub 

cal/ g 

-591.5  (c) 

Q;=+120 

K1 

2  2  n 
(kg/m-s)  /Pa 

1  K1=1.387xl0-4  (g) 

K°=9.619xl0~n  (g) 

nl 

- 

1.8  (f) 

1.8 

E1 

kcal/mol 

42  (f) 

• 

30. 

Tf 

K 

1982-0. 67T  (d) 

s 

1400 

Ts-RANGE* 

K 

835-910  (c) 

845-935 

#  All  actual 

calculations  were 

performed  using  the 

SI 

system. 

NOTES:  (a)  AMCP  706-177  Eng.  Design  Handbook  [35] 

(b)  Price,  Boggs  and  Derr  [33] 

(c)  Guirao  and  Williams  [32],  calculated  values 

(d)  See  comments  following  Table  2 

(e)  Taken  approximately  from  calculated  values  for  RDX  [31J 

(f)  Longwell  and  Wise  [36] 

(g)  See  Appendix  B  for  calculations 

(*)  2  MPa-10  MPa 
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The  following  comments  pertain  to  particular  values  of  parameters 
used  in  the  present  model  for  AP  as  shown  in  Table  1. 


from 


The  values  of  A  ,  E  in  the  pyrolysis  formula  were  calculated 
s  s 


JU  r  =  .AiCMpJ  - 


using  the  following  two  data  points  (r,  T  ): 

s 


p  (MPa) _ r  (cm/s)  *  Ts(K)»* 

2  0.3  835 

10  1.0  910 


*  Price,  Boggs  and  Derr  £33j,  r  vs  p  (measured) 

**  Guirao  and  Williams  [32j,  T  vs  p  (calculated) 

s 

6  2 

Thus, E  *  24.24  kcal/mol,  and  A  «=  1.29  x  10  g/cm  -s  are  obtained, 
using  a  sample  density  of  1.94  g/cc. 

The  value  of  Q*  s  444.6  cal/g  was  calculated  from  the  global 
enthalpy  balance. 


cc  (Jl  -To  )  ■*-  Cf.  (Tf  -Tt)  — Q*  -  o 

with  C  =  0.5,  C  *  0.3,  T  =  300K,  and  Tg  =  873K  (at  p  =  4.5  MPa, 
reference  point)?  Using  t$is  value  of  Q*,  one  may  obtain  the  linear 
dependence  for  Tf (T  )  shown  in  Table  1.  It  should  be  noted  that  the 
"surface  heat  release"  according  to  Eq.  (12)  is  thus 

Qg  —  (l— 5=1 

using  0*  h  =  -591.5  cal/g  (from  Guirao  and  V?illiams  C32])  and 
G  =  0.3fU  This  value  of  Q*  is  quite  close  to  that  inferred  by  BDP, 
namely, 120  cal/g. 
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TABLE  2 


HMX  INPUT  DATA 


PRQPEKTC 

UNITS# 

PRESENT 

MODEL 

BDF  MODEL 

PRES  EOT  CALCULATION 

A 

c 

cal/crr-s-K 

7xl0"4  (a) 

- 

C 

c 

cal/g-K 

0.43  (a) 

0.3 

pc 

g/cc  (pressed  pellets) 

1.88 

- 

A 

c 

Vs 

1019,7  (b) 

- 

E 

C 

kcal/mol 

52.7  (b) 

- 

n 

1  (b) 

__ 

c 

A 

s 

kg/m2-s 

2 . 09xl07  (c) 

5xl010 

Es 

kcal/mol 

19.59  (c) 

50. 

*g/cp 

kg/m-s 

5xl0~5  (c) 

10“3 

c 

p 

cal/g-K 

0.43=C  (c) 

c 

0.3=C 

c 

Qsub 

cal/g 

-141.6  (f) 

Q*=+225 

* 

Q1 

cal/g 

+453  (e) 

Q*=+667.5 

Tf 

K 

1353  (e) (g) 

3275 

T  -BANGE* 
s 

K 

650-767 

1085-1235 

E1 

kcal/mol 

46.2  (d) 

50=E 

s 

2  2  nl 

(kg/m -s)  /Pa  1 

1. 017xl0~7 

(h) 

4.34xKf12  (h) 

nl 

- 

1*4  (h) 

2 

#  Actual 

calculations  were  performed 

using  the 

SI 

system. 

NOTES : 

(a) 

AMCP  706-177  Eng.  Design  Handbook  [35] 

(b) 

Robertson  [37] 

(c) 

Ben-Reuven  and  Caveny  [38J 

(d) 

Shaw  and  Walker  [39] 

(e) 

Ben-Reuven  [31] 

* 

(f) 

Rosen  and  Dickinson  [40J 

(g) 

This  is  only  the  nearf ield-end  temperature, 
final  flame  temperature  (3100  K  at  20  atm) . 

(h) 

See  calculations  in  Appendix  B 

- 

(*) 

1  MPa-25.4  MPa 
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1 . 3  Discussion  of  Monopropellant  Simulation  Results 

The  formulation  developed  herein  was  used  along  with  the  BDP 
model  to  obtain  parallel  simulation  of  burning  rate  and  pressure 
sensitivity  data  for  both  AP  and  HMX. 

The  results  of  burning  rate  vs  pressure  are  shown  in  Fig.  5. 

The  experimental  data  were  taken  from  several  sources  J[  33 , 34  1  ir-ost 
of  which  were  unavailable  at  the  time  the  BDP  model  was  developed. 

For  this  reason,  both  models  were  adjusted  (single  parameter  fitted, 
pertaining  to  the  chemical  kinetics  constant)  at  a  single  experimental 
reference  print  (p,  m) ,  as  explained  in  Appendix  B.  The  plots  in 
Fig.  5  demonstrate  that  both  models  simulate  well  the  AP  burning  rate, 
while  the  present  model  is  clearly  better  for  HMX  in  the  pressure 
range  below  10  MPa  (100  atm.)  . 

The  same  HMX  data  set  was  previously  used  to  obtain  the 
parameters  (c,  B)  in  the  following  burning  rate  correlation, 


r 


(30) 


derived  in  Ref.  38,  and  reproduced  in  Fig.  6.  Evidently,  simulation 
by  Ec.  (30)  is  superior  to  both  BDP  and  the  present  model.  This  is 
not  surprising,  since  Eq.  (30)  incorporates  physical  details  pertinent 
to  the  complex  nature  of  nitramine  deflagration  (namely,  the  concept 
of  near-field  and  far-field,  each  dominated  by  a  distinct  chemical 
reaction  mechanism) .  These  details  are  completely  missing  from  the 
BDP  model,  and  were  omitted  from  the  present  model  (which  is  based 
entirely  on  the  near-field)  to  facilitate  the  single  delta- 
function  flame  approximation. 

The  pressure  sensitivities  due  to  each  model  are  compared  in 
Fig.  7,  where  the  trends  typical  of  AP  (n  decreasing  with  increasing 
p)  and  HMX  (n  increasing  with  increasing  p)  are  demonstrated.  Also 
shown  is  the  remarkable  capability  of  the  BDP  model  to  obtain  both 
increasing  and  decreasing  n (p)  with  the  same  functional  relationship, 
as  mentioned  earlier.  The  magic  of  this  achievement  is  entirely  due 

to  the  function  f  (T  ),  given  in  Eq.  (26). 

B  S 

Accordinq  to  Eqs.  (1)  and  (3),  the  BDP  model  provides  q  =  -lg(zR). 

From  Eq.  (3)  , 


(31) 


where  b.  =  C  /Q*  and  b,  =  (C  T  +  0*)/C>*.  Thus,  Eo.  (26)  yields 
0  .c  .f  1'cO  s'  f  1 

an  expression  in  terms  of  q. 


(32. a) 


20 


2 

where  the  multiplicative  constant  1/b^  has  been  omitted,  without 
loss  of  generality.  Note  that  is  positive  for  net  surface 
exothermici ty ,  as  assumed  by  BDP.  In  general,  0  <  q  <  1;  at  both 
ends  of  this  domain,  f^  — o  ,  and  the  single  extremal  point,  q*, 
is  defined  by 

—  0  (32. b) 

/s 

where  f^tq*)  is  a  minimum.  Now,  from  the  denominator  of  Eq.  (25. a) 
one  can  readily  see  that  whenever  €*  increases,  n  will  decrease, 
and  vice  versa  for  decreasing  f^.  Thus, the  question  of  n(p) 
progressivity  or  regressivity  in  the  BDP  framework  amounts  to  proper 
choice  of  parameters.  The  actual  parametric  ranges  employed  by  BDP 
are  as  follows. 

For  AP,  the  range  is 

T£(K)  P  (MPa)  q  ^  7 g 


830  1.58  0.186  4.50  .  1.68 

970  18.4  0.386  5.23  0.95 

for  bQ  =  1.43  x  10  3  1/K  and  b  =  1.0;  q* (b  )  =  0.2137  and  T*(q*)= 
849K,  while  p*  =  2.34  NPa.  At  this  point  n(p*)  is  maximal?  thus, 
n  will  increase  for  p  <  p*  and  decrease  for  p  >  p*,  as  shown  in  Fig. 
7.  This  behavior  is  entirely  extraneous  to  the  physics  built  into 
the  EDP  model,  and  is  not  collaborated  bv  experimental  observation. 

For  HN.X,  the  range  is 

Tg(K)  p(KPa)  q  Zg 

1070  0.7  9  x  10~3  5.46  4.71 

1235  21.6  8.3  x  10~2  1.49  2.49 

for  bp  =  4.49  x  10~4  1/K  and  b  =  0.472;  qMb^  =  0.1478  ,  an  order 
of  magnitude  higher  than  the  HFX  q-range.  Thus,  T*(q*)  =  1379  K, 
and  n(p)  will  be  progressive  throughout  the  range  of  p  considered, 
as  shown  in  Fig.  7. 
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Obviously,  variation  for  AP  is  moderate  in  comparison  with 
that  of  HMX ,  since  for  AP  the  q-region  straddles  the  f^-min  point. 

For  HMX,  the  q-region  is  assigned  (by  proper  choice  of  the  range 
and  other  parameters)  away  from  the  f^-min  point,  so  that  f^  variation 
is  more  pronounced.  Consequently,  one  may  expect  the  BDP  model  to 
obtain  steeper  n(p)  and  <$  (p)  variation  for  HMX  than  for  Ap  in  the 
same  pressure  range.  Thii  is  demonstrated  in  Figs.  7-11. 

Regarding  temperature  sensitivity,  both  the  BDP  model  and  the 
present  model  seem.  to  perform  rather  poorly  when  compared  to  the 
experimental  data  of  Price,  Boggs  and  DerrC33l.  The  associated  data 
reduction  for  AP  and  HMX  is  shown  in  Figs.  8  and  9, respectively . 
Calculated  and  experimental  values  of  c5  vs  p  are  depicted  in  Figs.  .10 
and  11.  For  AP,  Figs.  8  and  10  demonstrate  nonmonotonous  behavior 
(with  a  minimum  between  5  and  6  MPa),  but  relatively  small  overall 
variation:  100  d  between  0.2  and  0.3  1/K,  for  1  <  p  <  10  MPa.  In 

the  same  pressure^range,  d  for  HMX  is  monotonously  decreasing,  with 
a  relatively  large  overall  Variation:  100  <6  varies  between  +0.4  and 
-0.04  1/K,  as  shown  in  Fiqs.  9  and  11;  crossing  to  negative  values 
occurs  roughly  at  7  MPa.  Note  that  the  BDP  results  for  HMX,  although 
following  the  general  trends  at  low  pressures  in  Fig.  11,  are  still 
an  order  of  magnitude  too  high. 

The  folloving  conclusions  can  be  drawn.  Both  the  BDP  model  and 
the  present  model  perform  reasonably  well  in  simulation  of  burning 
rate  data,  with  a  single  parameter  adjustment.  Both  models  obtain 
explicit  solutions,  although  the  present  model  is  somewhat  mere 
involved  alaebraically ,  on  account  of  the  input  reactant  fraction, 

G(T  ,  m?  T^) ,  which  is  completely  absent  from  the  BDP  formulation. 

BotFi  models  perform  rather  poorly  in  simulation  of  pressure-  and 
temperature-sensitivity  of  burning  rate. 

The  major  difference  between  the  models  is  quite  apparent.  The 
present  model  can  be  derived  rationally  from  basic  principles  of 
combustion  theory  (using  both  energy  ana  species  conservation 
equations)  and  involves  functional  details  (such  as  subsurface 
decomposition  process  and  chemical  kinetics  data)  specific  to  the 
monopropellant  under  consideration.  In  contrast,  the  BDP  model 
incorporates  an  intuitive  argument  (extraneous  introduction  of 
chemical  kinetics  effects) ,  while  excluding  specific  physical 
details;  consequently,  the  BDP  formulation  is  not  in  agreement  with 
known  combustion  theory,  and  characteristics  of  each  monopropellant 
must  be  entered  by  proper  choice  of  parameters  and  T  -range.  This 
parameter  selection  is  often  tenuous,  as  evident  from  the  high  values 
of  T  imposed  for  HMX  (which  seem  more  adequate  for  the  outer  end 
of  tfie  gaseous  nearfield,  where  primary  decomposition  is  complete)  , 
or  from  the  n(AP)-max  obtained  inadvertently. 

The  delta-function  approach  developed  herein  is  admittedly 
oversimplified  (for  HMX  in  particular),  and  seems  too  crude  for 
anything  but  burning  rate  vs  pressure  simulation.  It  shows, 
nevertheless,  that  a  reliable  alternative  to  the  BDP  monopropellant 
model  can  be  derived,  and  has  potential  for  impr overmen t,e.g.  , better 
subsurface  decomposition  formulation. 
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1 . 4  The  BDP  Composite  Propellant  Model 


This  section  addresses  the  problems  associated  with  the  BDP 
combustion  model,  specifically  pertaining  to  composite  propellants, 
as  shown  in  the  three  parts  of  Fig.  3. 

The  BDP  model  employs  two  consecutive  flames  in  its  physical 
picture,  cf  Fig.  3. a,  namely,  the  "monopropellant  flame"  and  the 
"final  diffusion  flame."  Both  are  assumed  adiabatic,  and  their 
respective  heat  feedbacks  to  the  surface  are  simply  added;  in  other 
words,  the  two  flames  are  assumed  independent  of  each  other,  or 
decoupled.  Obviously,  conductive  heat  transfer  from  the  external 
(diffusion)  flame  must  pass  through  the  premixed  monopropellant  flame 
region.  An  approximate  double  delta-function  flarre  analysis  (in  which 
the  external  flame  position  is  assumed  to  be  far  downstream  from  the 
inner  flame)  indicates  that  appreciable  nonlinear  coupling  between 
the  external  and  inner  flames  should  prevail.  The  following  results 
demonstrate  this  point. 

Suppose  that  the  external  flame  has  an  exothermicity  Q*, 
comparable  to  that  of  the  inner  flame,  and  is  placed  at  a 
dimensionless  distance  of  z^  from  the  surface,  such  that  z2»z^. 

Thus,  one  may  consider  the  external  flame  as  a  small  perturbation 
in  the  region  close  to  the  surface.  The  thermal  enthalpy  at  the 
position  of  the  inner  flame  is 

01  =  ©fed  =  A+ <33> 

where  the  last  term  in  parantheses  represents  the  dimensionless  heat 
feedback  to  the  surface,  and 

^z~C  qv)  ^  0.z Cf^Tz.)/ to2,  (34) 

Now,  relative  to  the  unperturbed,  single  delta-function  inner  flame, 
the  perturbation  is 

o  Se,  =ei-e?=cz£*zC<£*‘-\ )  „S) 

where  (  )  denotes  unperturbed  properties;  note  that  z^  =  has  been 
assumed  (as  inferred  from  detailed  calculations,  cf  Fig.  6  of  Ref. 

38),  as  well  as  A  =  A  .  This  leads  to  a  perturbation  of  the  inner 
flame  temperature,  T^(z^) 

-  CQ*/qV)  sa, 
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and  consequently,  through  the  eigenvalue  relationship  in  Eqs,  (15)  , 
(16) ,  the  associated  burning  rate  perturbation  is 


\wi/yr\°  a/  Srj/T/ 


(36) 


with  /2T  )  0(10),  demonstrating  appreciable  amplification. 

This  nonlinear  coupling  effect  is  completely  absent  from  the  BDP 
composite  propellant  analysis,  since  adiabatic  values  of  =  T^ 
are  specified,  and  the  surface  enthalpy  balance  (where  all  heat 
feedback  components  are  linearly  superimposed)  is  used  to  calculate 
the  overall  burning  rate.  It  is  therefore  indicated  that  the  BDP 
model  is  deficient  in  this  respect  of  external/inner  flame  coupling, 
even  if  the  physical  model  of  Fig.  3.b  is  accepted. 


Another  puzzling  aspect  of  the  BDP  composite  model  is  that  the 
parabolic  (or  conical)  diffusion  flame  in  the  physical  picture,  cf 
Fig.  3,  is  transformed  into  a  planar  sheet  parallel  to  the  propellant 
surface  in  the  physical  model,  Fig.  3.b.  Then,  in  the  analytical 
model,  BDP  assigns  a  value  of  exp(-c*zD)  to  the  heat  feedback  due 
to  this  flame,  where  z  is  the  dimensionless  flame  height, 
proportional  to  the  actual  height,  X*^,  and  c  denotes  a 
proportionality  constant . 


In  an  attempt  to  see  whether  such  mean  flame  height  is 
analytically  plausible,  an  approximate  calculation  of  the  mean  heat 
feedback  is  carried  out,  for  a  conical  flame  with  a  height  of  z^. 

The  configuration  is  shown  in  Fig.  12.  The  prevailing  assumption 
is  that  the  heat  feedback  contribution,  due  to  each  circular  element 
of  the  conical  flame  surface,  is  proportional  to  exp(~z),  where  z 
denotes  the  element  height.  The  contributions  from  the  entire  cone 
surface  are  then  integrated,  and  the  corresponding  mean  heat  feedback 
to  the  oxidizer  crystal  surface  is 


— 

nez 


,  Jes2Tr(ec-2ty>Otno<<fe 
c  0 

=  2  Op- H- 


(37) 


where  R  is  the  crystal  radius.  The  expression  obtained  herein  for 
the  mean  heat  feedback ,  -q  (zD) ,  clearly  precludes  simple  reduction 
to  the  form  used  by  BDP,  where  z appears  only  as  the  exponential 
argument. 
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The  f*inal  point  to  be  raised  in  this  conjunction  concerns  the 
validity  of  the  so-called  final  diffusion  flame.  EDP  postulated  that 
this  flame  is  diffusion  controlled  under  all  conditions,  and  without 
reference  to  particular  chemical  kinetics  length  or  time  scales. 

This  means  that  chemical  kinetics  relaxation  times  are  typically 
much  shorter  than  diffusion  times.  Clearly,  the  reactant  streams 
(emanating  separately  from  binder  pyrolysis  and  from  oxidizer  primary 
decomposition)  are  non-premixed  or  initially  separate?  further,  the 
expected  temperature  range  should  be  quite  high,  close  to  the 
adiabatic  flame  temperature  of  the  oxidizer/binder  configuration. 

All  this  seems  to  favor  fast  reaction  kinetics  and  diffusion- 
controlled  sheet  flame.  However,  calculations  carried  out  with 
chemical  reaction  length  scales  relevant  to  nitramine  combustion 
demonstrate  that  some  typical  secondary  reactions  have  length  scales 
in  the  order  of  100  microns  and  more.  This  is  shown  in  Fig.  13;  these 
scales  are  similar  to  typical  hydrocarbon-oxidation  reactions 
expected  between  binder  chain  fragments  and  oxidizer  primary 
decomposition  products  for  both  AP  and  HMX.  Moreover,  oxidation 
reactions  involving  NO,  a  typical  HMX  decomposi tion  product, are 
expected  to  be  even  slower.  This  strongly  indicates  that  in  some 
cases,  specific  to  both  propellant  configuration  and  pressure  regime, 
the  final  flame  might  be  chemical-kinetics  controlled,  or  that 
kinetics  and  diffusive  processes  might  have  comparable  length  and 
time  scales. 
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1 .5  Conclusions  and  Recommendations 

The  foregoing  results  point  to  several  deficiencies  in  the  BDP 
combustion  model  for  composite  propellants.  These  are:  (1)  The 
monopropellant  (premixed)  flame  model  does  not  agree  with  known 
combustion  theory.  (2)  The  analytical  composite  model  combines 
linearly  heat  feedback  from  three  different  flame  components  without 
any  consideration  of  thermal  coupling,  a  prevalent  nonlinear  effect. 
(3)  Heat  feedback  from  the  parabolical  or  conical  flame  surface  is 
not  adequately  represented.  (4)  The  final  flame  is  assumed 
diffusion-controlled  under  all  conditions,  regardless  of  relevant 
chemical  kinetics  or  pressure  regime.  Some  of  these  defects  could 
have  been  readily  corrected.  It  can  only  be  lamented  that  such  a 
large  amount  of  mental  energy  went,  over  a  period  of  ten  years,  into 
considerations  of  granularity,  modality,  concave/convex  particle 
surfaces,  distinct  vs  uniform  binder/oxidizer  surface  temperatures, 
etc.,  —  all  in  support  of  a  semi-empirical  correlation  (exponential, 
not  polynomial)  which  does  not  live  up  to  its  own  physical  model. 

Evidently,  this  analysis  indicates  two  main  routes  for  future 
solid  propellant  combustion  research.  One,  with  an  immediate 
engineering  utility,  would  be  to  use  directly  Chebichev  polynomials 
of  the  first  kind  to  simulate  burning  rate  vs  pressure,  with 
adjustable  coefficients  to  account  for  granularity,  modality,  etc. 
These  must  work,  since  they  involve  exp(-z): 


-to O 


w 


Chebichev  polynomials ,  order  n=l  through 
5  plotted  against  normalized  independent 
variable,  w. 
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where  Jn  denote  the  associated  Bessel  functions  of  the  first  kind, 
arid  w  is  a  normalized  independent  variable*  These  polynomials  are 
routinely  available  in  any  respectable  minicomputer  software  package, 
whereby  a  large  number  of  coefficients  can  be  adjusted  with  minimal 
effort. 

The  second  way,  by  which  physical  insights  might  be  obtained, 
would  be  through  detailed  analysis,  using  explicit  chemical  mechanisms 
and  kinetics  data  in  a  comprehensive  model,  which  may  probably  require 
numerical  solutions  of  differential  systems. 
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Fig.  1  The  physical  concept  of  burning  composite 
propellant,  due  to  BDP. 
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Fig.  2  The  Combustion  Fraction,  an  indicator  of  the  relative  importance 
of  combustion  processes  to  various  propellant  burning  models.  The 
trend,  obtained  by  least  squares  fit  (excluding  the  model  by  King  [ 18 ] ) 
is  obviously  negative. 
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Fig.  3  Evolution  of  the  BDP  analytical  model  from  the  physical  model 
and  the  physical  picture  imposed  by  the  authors.  One  major  question 
in  the  critique  herein  is  whether  the  three  models  are  self-consistent . 
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Fig.  4  The  delta  function  approximation  for  the  gas  phase,  used  to 
develop  the  present  simplified  model.  The  delta  function  is  imposed  only 
in  the  dimensionless  physical  coordinate  system,  (0,B) . 
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Fig.  5  Comparison  of  burning  rate  vs  pressure  for  HMX  and  AP  by  both 
BDP  and  the  present  model  (denoted  MBR) ,  showing  the  latter  to  obtain 
better  correlation  for  HMX  below  10  MPa. 
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Fig.  6  Nitramine  burning  rate  vs  pressure  -  the  same  HMX  data  used 
in  Fig.  5  is  correlated  by  Eq.  (30)  ,  demonstrating  superiority  over 
both  BDP  and  the  present  model  calculations. 
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Fig.  7  Burning  rate  pressure  sensitivity,  n,  calculated  for  both  AP 
and  HMX.  Both  BDP  and  MBR  models  yield  poor  results  for  HMX,  which  is 
correlated  by  use  of  Eq.  (30)  . 
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Fig.  8  Regression  rate  vs  initial  temperature  for  AP,  with  pressure  as 
parameter.  Correlation  of  the  experimental  data  by  Price,  Boggs  et  al 
[33]. 
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Regression  rate  vs  initial  temperature  for  HMX,  with  pressure  as 
parameter.  Correlation  of  the  experimental  data  by  Price,  Boggs 
et  al  [33 j  .  ^ 
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Fig.  10  Calculated  and  experimental  temperature  sensitivity  vs  pressure, 
for  AP.  Experimental  points  are  summarized  from  Fig.  8. 
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100  x  OL  (l/K) 


Fig.  11  Same  as  Fig.  10,  for  HMX.  The  experimental  trend  is  followed 
by  the  BDP  model  to  some  extent,  but  an  order  of  magnitude  higher. 
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Fig.  12  The  conical  flame-sheet  approximation,  to  calculate  an  integrated 
mean  heat  feedback  to  the  surface. 
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PRESSURE,  MPa 


Fig.  13  Comparison  of  chemical  reaction  and 
diffusive  length  scales,  relevant  to  AP  and  RMX . 

Even  for  the  oxidizer  (fuel  flame  of  AP,  calculated 
by  Steinz  f6])^  the  reactive  and  diffusive  scales  may 
be  comparable,  so  the  "final  diffusion  flame"  concept 
of  BDP  should  be  re-examined. 
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II.  A  COMPREHENSIVE  ANALYTICAL 
MODEL  FOR  NITRAMINE  DEFLAGRATION 


1 . 1  Background 

The  foregoing  .section  indicates  that  the  complexity  of  nitramine 
monopropellant  deflagration  has  not  been  considered  in  many  composite 
propellant  modeling  efforts.  Evidence  to  this  complexity  can  be  found 
in  abundance  in  the  nitramine  chemical  decomposition  literature,  a 
sample  of  which  is  given  in  Refs.  1-22;  these  works  have  been  thoroughly 
reviewed  in  Ref.  23  and  from  the  basis  of  a  detailed  analytical  model 
[2k,  25];  the  present  work  is  a  direct  extension  of  this  model.  Further 
reviews  of  published  chemical  analyses  were  performed  by  Schroeder  [26-28]. 

The  main  objective  of  the  present  analysis  is  to  enable  burning  rate 
prediction,  i.e.,  calculation  of  m  (p;Tq)  over  a  wide  pressure  range. 

Since  the  near-field  in  the  gas  phase  has  been  treated  in  some  detail 
previously  [23-25],  in  terms  of  two  global  reactions,  a  natural  extension 
would  be  to  incorporate  the  entire  gaseous  flame  (including  the  far-field) 
and  to  improve  the  melt  phase  model.  This  entails  a  larger  number  of 
secondary  reactions  and  chemical  species  than  previously  considered.  Fur¬ 
ther,  in  order  to  calculate  the  burning  rate  independently,  full  closure 
of  the  gas/melt  interface  conditions  is  required;  this  is  provided  by 
incorporation  of  a  nonequilibruim  evaporation  law.  This  part  of  the  study 
is  motivated  by  the  desire  for  better  understanding  of  the  nitramine 
combustion  mechanism,  and  ways  to  affect  changes  in  the  burning  rate, 
ultimately.  The  capability  to  compute  the  flame  field  structure  (species 
and  temperature  profiles)  as  well  as  temperature  and  pressure  sensitivities 
of  the  burning  rate  would  enable  comparison  with  experimental  data  for 
verification. 

This  leads,  of  course,  to  greater  computational  complexity  and  points 
up  the  relative  lack  of  reliable  chemical  kinetics  data.  In  this  respect, 
the  analytical  model  developed  herein  could  serve  to  test  various  relevant 
chemical  mechanisms  and  kinetics  data,  despite  the  higher  degree  of  un¬ 
certainty  involved. 

The  following  sections  describe  the  elements  in  the  analysis  which 
are  new  relative  to  Refs.  23-25:  An  extended  (highly  idealized)  chemical 
mechansim  for  the  gas  phase,  the  melt  phase  model  including  decomposition 
gas  bubbles,  the  gas  phase  model  extended  to  include  the  far-field,  and 
the  nonequilibrium  evaporation  law. 
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1.2  CHEMICAL  MECHANISM  AND  KINETICS 


Nitramine  deflagration  probably  involves  two  major  sets  of 
chemical  reactions:  (1)  primary  decomposition,  first  order  overall, 
(in  both  condensed  and  gas  phases)  in  which  the  molecular  substance 
breaks  up  to  form  intermediate  reactants,  which  is  relatively  fast 
and  would  occur  close  to  the  propellant  surface,  and  (2)  secondary 
reactions  between  the  intermediate  species,  typically  second  order 
overall,  relatively  slow  (compared  with  primary  decomposition) ,  by 
which  the  final  composition  and  the  final  temperature  evolve.  This 
conceptual  division  into  two  rather  diverse  families  of  reactions 
is  based  on  reported  observations  of  the  mechanism,  products  and 
kinetics  of  nitramine  decomposition£l-223  ;it  tends  to  be  strongly 
supported  by  the  measured  shape  of  burning  rate  vs  pressure  curve, 
which  exhibits  a  progressive  pressure  exponent  n(p).  For  detailed 
discussions  of  the  mechanism  and  kinetics  of  nitramine  decomposition, 
the  reader  is  referred  to  Refs.  23,24  and  25  . 


The  chemical  reactions  considered  in  the  present  analysis  are 


as  follows: 

RDX(L)  -*-3CH20+9/4N20+3/4N0+3/8N2,  Q=51.4 

kcal 

(R.l) 

RDX  (G)  -v  N02+N20+3CH2CH-3/2N2 , 

Q*=106.8 

kcal 

(R.2) 

5/7CH20+N02-^N0+3/7CX)+2/7C02+5/7H20,  Q=46 

kcal 

(R.3) 

ch2o+n2o  -*-h2o+co+n2  , 

Q  *  76.1 

kcal 

(R.4) 

ch2o+no  — *-h2o+co+i/2n2, 

Q  =  78.1 

kcal 

(R.5) 

CO  +  h2o-*>co2  +  h2. 

Q  =  9.9 

kcal 

(R.6) 

co2  +  H2_*‘co  +  H20' 

Q  =  -9.9 

kcal 

(R.7) 

where  Q  ">  0  denotes  exothermicity.  All 

the  reactions  (R. 

2)  through 

(R.7)  can  occur  simultaneously  at  any  point  in  the  gaseous  flame 
field.  The  global  mechanisms  postulated  for  reactions  (R.4)  and  (R.5) 
are  purely  conjectural,  and  merely  intended  to  produce  observed  final 
products.  The  relevant  chemical  kinetics  data  is  given  in  Table  1. 

It  should  be  emphasized  that  the  list  of  chemical  reactions  herein 
is  in  no  way  conclusive.  Many  more  reactions  may  be  relevant  to  the 
gas  phase  during  nitramine  deflagration;  further,  some  intermediate 
products  might  have  a  strong  influence  upon  the  far-field  structure, 
the  final  composition,  and  final  temperature  attained.  The  algorithm 
constructed  in  this  study  for  solution  of  the  conservation  equations 
in  the  gas  phase  anticipates  progress  in  the  detailed  knowledge  of 
relevant  mechanisms  and  chemical  kinetics  data.  Incorporation  of 
additional  chemical  reactions  in  the  model  involves  merely 
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modifications  of  the  input  data  set,  as  the  reaction  terms  are  handled 
in  a  completely  generalized  manner  in  the  algorithm.  For 
accomplishment  of  this  flexibility,  the  list  of  chemical  species  (and 
the  specific  thermophysical  data  associated  with  them,  such  as 
specific  heats,  conductivities  and  viscosities,  molecular  weights, 
etc) ,  is  more  extensive  than  would  be  required  by  the  foregoing  list 
of  reactions,  and  is  as  follows: 

nitramine  (RDX  or  HMX)  ,  H20,  CH20,  CC>2,  CO,  N2, 

NO,  N20,  N02,  HCN ,  HN02,  H2,  H,  HO,  and  (-CH^) 

Evidently,  reactions  by  which  HCN  +  HONO,  or  HONO  4-center 
elimination  occur  were  excluded  from  the  present  analysis. 
HONO-elimination  was  suggested  by  Shaw  and  Walker [20]  as  a  possible 
initial  step  of  nitramine  decomposition,  later  followed  by  Schroeder, 
£26-“ 26]  and  recently  by  McGuire  and  Traver  [29J ;  all  these  works  involve 
theories,  not  measurements.  Experimentally,  the  Electron  Spin 
Resonance  (ESR)  spectroscopy  conducted  by  Beyer  and  Morgan  [3oJ 
indicated  the  presence  of  HCHN  and  N02  during  decomposition  of  RDX 
and  HMX;  these  may  lead  to  subsequent  HCN  +  N0„  reaction;  however, 
HONO-elimination  is  not  a  necessary  step  and  HCHO  production  is  still 
present.  It  must  be  stressed  that  identification  of  large  molecular 
fragments  of  74  and  148  amu  by  Goshgarian  £l6]  and  by  Farber  and 
Srivastava, £31 J  by  mass  spectroscopy  of  decomposing  RDX  and  HMX  do 
not  lend  any  evidence  to  the  HONO-elimination  schemes;  rather,  further 
support  is  obtained  for  the  widely  observed  HCHO,  NO  and  NO 
decomposition  products.  For  these  reasons,  the  reactants  HCN  and 
HONO  are  incorporated  in  the  data  set  of  the  present  analysis,  but 
the  corresponding  reactions  were  suspended  until  further  experimental 
evidence  points  out  their  relevance.  As  to  HCN  oxidation  by  N02, 
the  high-temperature  and  pressure  (but  dilute  mixture)  measurements 
by  Fifer  and  Holmes  £32] indicate  that  this  overall  reaction  is  slower 
by  1  to  3  orders  of  magnitude  than  the  parallel  N02  +  HCHO  reaction 
in  the  range  1000-1500k  typical  to  the  end  of  the  gaseous  near  field 
in  the  deflagration  process,  which  means  that  the  HCN  +  NO  oxidation 
does  not  influence  the  heat  feedback  to  the  propellant  surface 
appreciably. 


1.3  ANALYSIS 


Generally,  the  physical  model  of  nitramine  monopropellant  deflagration 
herein  is  quite  similar  to  the  models  discussed  in  detail  by  Ben-Reuven  et  al; 
(23-25J  therefore,  only  those  elements  of  the  present  analysis  which  are  distinct 
will  be  presented.  In  particular,  (1)  incorporation  of  gas  bubbles  in  the  melt 
layer,  (2)  extension  of  the  analysis  to  include  the  far-field  in  the  gas  phase, 
and  (3)  addition  of  a  nonequilibrium  nitramine  evaporation  law,  to  facilitate 
independent  calculation  of  the  mass  burning  rate,  m. 

The  physical  models  for  the  gaseous  deflagration  wave  and  the  melt  layer 
are  depicted  in  Figs.  1  and  2, respectively.  The  analytical  method  employed 
is  to  define  and  solve  separately  well— posed  problems  in  each  region  (melt, 
and  gas  phase),  and  then  match  the  solutions  at  the  melt/gas  interface  by 
satisfying  the  available  mass,  species  and  energy  conservation  constraints. 
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1.4  CHEMICALLY  REACTING  MELT  LAYER  WITH  GAS  BUBBLES 


The  presence  of  gas  bubbles  in  the  melt  layer  of  deflagrating 
neat  nitramines  (RDX  and  HMX)  has  been  inferred  from  electron 
microscope  photography  of  extinguished  sample  surfaces  [KQ • Their 
relevance  to  the  actual  deflagration  process  was  further  indicated 
in  a  previous  theoretical  analysis!^ where  the  effect  of  bubbles  was 
identified  a  posteriori  by  an  approximation  (which  allowed  a  unique 
correlation  between  the  imposed  global  heat  of  reaction  in  the  layer 
and  pressure) .  The  melt  layer  analysis  herein  incorporates  gas 
bubbles  a  priori,  through  a  two-phase  flow  treatment. 

It  should  be  emphasized  that  even  with  1%  subsurface 
decomposition  in  which  bubbles  are  formed,  the  local  mean  density 
would  be  appreciably  lower  than  the  neat  liquid  density;  of  course, 
this  follows  from  the  great  diversity  in  intrinsic  densities  at 
moderate  pressures:  <^*  Similarly,  since  thermal  conductivities 

are  related  approximately  as  0(10),  the  porosity,  due  to 

the  presence  of  bubbles,  can  effectively  decrease  the  conductive  heat 
transfer  in  the  layer.  In  comparison,  the  available  data  show  that 
Cc  =  0^  (gas)  in  the  range  of  temperatures  considered. 

The  following  simplifying  assumptions  have  been  employed: 

(1)  Bubbles  are  small  relative  to  the  overall  layer  thickness. 

(2)  Bubbles  are  small  enough  such  that  thermal  relaxation  within  the 

bubbles  and  between  the  bubbles  and  the  liquid  can  be  considered 
instantaneous,  and  a  uniform  local  temperature,  T^  =  T  prevails  at 
each  point.  (3)  Mean  intrinsic  thermophysical  properties 
(conductivities  and  specific  heats)  are  uniform  throughout  the  layer, 
as  well  as  the  liquid  density,  <p  *.  (4)  Further,  the  pressure 

throughout  the  layer  is  uniform:  cdp/dy  «  0,  and  equal  to  the  externa 
pressure  prevailing  in  the  gas  phase  domain.  (5)  Although  thermal 
diffusion  is  incorporated,  chemical  species  diffusion  (due  to 
concentration  gradients)  is  assumed  negligible  in  the  layer;  note, 
however,  that  the  gas  may  percolate  through  the  liquid,  as  a  nonzero 
velocity-slip  is  permissible. 

For  two  interpenetrating  media,  the  total  volume  is  V  =  V + 

Vc.  The  porosity  or  M void-fraction”  is  then'\|/'=  V  /V,  while  tne 
overall  mixture  density  is  defined 


(1) 


where  £>*  =  const,  and  =  p  W  /R^T  are  intrinsic  liquid  and 
gas  densities,  and  Y  denotes  th§  overall  liquid  mass  fraction. 
According  to  the  foregoing  definition  of  porosity, 


(2) 
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Similarly, 


i  =  Vc/V  =  Vf6  /?* 

both  of  which  will  become  useful  later.  The  mass  conservation 
equations  arise  from  a  continuum  derivation  for  two  phase  flow,  such 
as  given  by  Nigmatulin  [Ml ; at  steady  state, 

-  3)  [0 -V)?6U6]=WC  ,3a, 

ac]=-ojc  ,3b) 

3 

Cx)  (kg/m  -s)  denotes  the  mass  exchange  rate  between  the  phases,  due 
to  liquid  nitr amine  decomposition  by  which  a  gas  mixture  evolves, 

^  =  A«.e^Tftl)(?6v)  ,4, 

Summation  over  the  two  mass  conservation  equations  yields  the  overall 
constraint:  dm/dy  =  d/dy(m^  +  m^)  *  0,  from  which 

0-v)?6«6tVf6uc=w(i1)=?tV=^1t  (5) 

determined  at  the  solid/liquid  boundary,  where  Y(-  y  )  =1,  and 
u  {-  y  )  *  r.  The  formal  similarity  between  the  foregoing 
differential  equations  for  liquid  and  gas  can  be  utilized,  so  that 
only  one  need  be  solved  (e.g.,  the  liquid  phase),  as  will  be  shown 
later. 


The  liquid  phase  'momentum  equation  in  the  layer  is 

jctj  [C'-xl/')?c't'Uci]+  (I =  -WccOc  (6, 

which  reduces  to  duc/dy  =  0 ,  or 

U-cOp  «  r  =  Cc7K.st  (7) 

where  assumption  (4)  and  Eq.(3.b)  have  been  applied.  No  useful 
information  is  obtained  from  the  gas  phase  momentum  equation,  since 
with  uc  *  r  and  *  fully  determined  by  T(y),  one  may  obtain 
algebraically  from  Eq.  (5) , 

wA£(TC^y)  <3> 

which  is  similar  to  the  definition  of  velocity  in  the  gaseous  flame 
field,  where  m  =  const  and  p  =  const  are  likewise  assumed. 
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Utilizing  assumptions  (2)  and  (3) ,  the  combined  energy  equation 
(Over  both  gas  and  liquid  phases)  is  written  in  terms  of  the  overall 
thermal  enthalpy,  hc  =  C0  (T  -  T^) : 


(9) 


where  =  C  =  C  =  const,  and 
B  c  p 

Ab  =  +  0  -1Jr)  ,10) 

the  overall  conductivity  is  now  obviously  variable  through  the  layer, 
through  its  dependence  upon  the  porosity.  It  should  be  emphasized 
that  the  layer  thickness  is  not  known  beforehand;  however,  three 
boundary  conditions  are  available  for  the  energy  equation,  viz., 

TC-Mw)  =  Tn  /  T(p)=rs 


lr,rCcC't«'To)  +  0»  3.  (ID 


which  help  overcome  this  difficulty.  The  following  coordinate 
transform  is  used: 


dy' 

As/mCe 


(12) 


This  transform,  and  the  dimensionless  variables 
and  X  s  re<^uce  differential  system 


TTs  (T-T  )/(T  -T  )  , 

r:  a  m  S.  m 

for  0  <  z  <  z  to: 
s 


(13) 


=  _<s 

35  • 


(14) 


-CIO')  =0  dZ/dZ(0)  =  [QZl-&(Tm-To)l/CcCTi-To), 

X(0)  =  1  .  (1 


where  §c  =  Q*/C  (T  ~Tm) »  00,3  (  and  zg  defined  by  Eq.  (12), 

with  upper  limit  or  0.  A  final  transform  So  *CT  as  independent 
variable,  with  <£>  =  •£-  d'C'/dz  and  z  as  dependent  variables, obtains, 
for  0  <r  <  1: 

djpldx  =  Sc  Wc  /(x-  <t> ) 
ds/dr  = 

<£>(0)  =  -  decided)  ?  z(0)  =  O 


(16) 

(17) 

(18) 
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while  X(  £;)  is  given  at  each  point  by  the  coupling  coefficient 

c*  =  4>&)/Qc  +  XCc)  =  ^  Hco)  +•  1  =  a~sr(19 , 

where  the  self-similarity  between  *  and  X  in  Eqs.  (13)  and  (14)  has 
been  exploited.  Of  course,  an  additional  integration  is  necessary 
to  determine  y  from  the  dimensionless  z. 

The  formulation  in  this  section  therefore  obtains  a  closure, 

by  which  all  of  the  layer  properties  of  interest  (including  Y,  dT/dy, 

§  0t  u  and  T|^)  can  be  caluclated  as  functions  of  temperature,  with 

(m,  Tg;  p,  Tq)  as  parameters.  The  major  physical  difference  between 

the  present  melt  layer  model  (with  bubbles)  and  the  previous  ones 

[2i  ,2S3  , (without  bubbles)  is  that  pressure  effects  are  incorporated 

explicitly  herein,  whereas  in  the  previous  models  solutions  could 

be  generated  for  specified  (m,  T  ;  T  )  data. 

s  u 

The  differential  system  of  Eqs.  (16)  through  (18)  was  integrated 
numerically,  using  a  Runge-Kutta-4  method. 
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1 . 5  THE  GAS  PHASE  REGION,  INCLUDING  THE  FAR -FIELD 


In  the  case  of  nitramine  deflagration  at  moderate  pressures, 
the  diversity  of  chemical  kinetics  length  scales  (between  rapid 
primary  decomposition  and  much  slower  secondary  reactions)  implies 
that  the  gas  phase  region  is  divided  into  two  distinct  parts:  the 
near-field,  adjacent  to  the  surface,  where  rapid  change  and  high 
gradients  occur,  and  the  far-field,  where  changes  in  composition  and 
temperature  are  much  more  moderate  in  comparison. 


Employing  standard  assumptions  (for  a  discussion,  see  Refs.  Z>S 
and  23),  the  dimensionless  conservation  equations  for  chemical  species 
mass  fractions  and  thermal  enthalpy  are,  in  the  region  0  <L^<00: 


AA 


dX/Sk,  - 


-6=1 


(20) 


J  =(j2-  fO 


dvi/ck^  -  d^ldf 


(21) 


3 

where  oO.  (mol/m  -s)  denote  overall  reaction  rates,  and  the 
dif f usivi^ies  are  9v/C  =  ^D,  equal  for  all  species,  but  allowed 
to  vary  with  both  loca?  composition  and  temperature.  The 
dimensionless  coordinate  £  and  the  thermal  enthalpy  are, 
respectively, 


$ 


[W 

o 


(22) 


h=  gv  fcftrHr' 

-TO 


(23) 


where  T°  is  the  standard  reference  temperature  and  Q*  denotes  the 
heat  of  primery  nitramine  decomposition  in  the  gas  phase  (J/kg) . 
Note  that  m(^)  =  const.  The  boundary  data  specified  are: 


Yj  (o+)  =  wo*)  -  «rs  >  -  YT )  <24a> 

=-  W-j  ,  Tp  5  dVj /^C°°)=dK/^C<w) =0.  (24b) 

Since,  as  mentioned  earlier,  the  gaseous  flame  field  involves 
diverse  length  scales,  the  following  transform  was  used,  for  practical 
purposes: 
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with  c>0  being  an  adjustable  constant,  which  allows  controlled 
streching  of  the  near-field  coordinate,  while  compressing  the  far- 
field  length  scale.  It  also  affords  sufficient  spatial  resolution 
in  the  far-field  (e.g.,  for  specified  £  =  1000  and  c  =  0.3, 

n  -  6.91  and  yi  ,  =  -1.2)  .  max 

[  max  Lmm 


Using  the  foregoing  transform,  the  gas  phase  conservation 
equations  can  be  expressed  as 


cLd> 

drr\ 


(26) 


yr\CU)t 


where  4?(7|)  represents  dimensionless  gas  phase  variables  (species 
mass  fractions  thermal  enthalpy)  and  R  denotes  the  sum  of  relevant 
chemical  reaction  terms,  as  in  Eqs.  (20),  (21).  For  T|  <  0  (the 
near-field) ,  the  formulation  gives  approximately 


ddplchr\  -  dYp/drf-  =  6  £.(40 


where  0<£<.  1  (but  where  reaction  terms  are  relatively  large), 
while  for  7\  p-  0  (the  far-field) ,  the  diffusion  term  is  relatively 
suppressed,  and  a  convective-reactive  flow  prevails: 

d^/ck\  ~  e7  £(.$)  (28) 

Note  that  for  large  values  of  y\,  the  chemical  reaction  rates  are 
expected  to  be  vanishingly  small,  compatible  with  the  boundedness 
constraint  imposed  in  Eq.  (24. b).  Independent  solution  manifolds 
can  now  be  generated  for  the  gas  phase  conservation  equations 
(transformed)  {2^3 , for  specified  data  (m,p)  and  4>(  "HnnJ  =  ^  * 

The  latter  data,  of  course,  are  not  known  a  priori.  Therefore,  the 
relevant  interface  conservation  conditions  at  the  melt/gas  interface 
must  be  satisfied  to  render  the  overall  solution  (carried  out 
separately  for  the  melt  layer  and  the  gas  phase)  physically 
meaningful.  For  chemical  species  and  thermal  enthalpy,  following 
Scala  and  Sutton  [363, 

<k (o+)  -  dcp/dk, (0+)  b<H0")  ( 29 ) 

where,  for  chemical  species 

4*  (.O')  =  Yj (0"),  =  lj>2. /vi  oo) 

and  for  thermal  enthalpy, 


Jp  (0-)  =  [  mce,  Cn  -T»-)  -  ^  (o') 

on 
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where  the  right-hand  sides  and  the  left-hand  sides  of  Egs.  (29)  are  fully 
determined  by  the  melt  layer  and  the  gas  phase  solutions,  respectively.  Those 
solution  manifolds  which  satisfy  the  interface  conservation  conditions,  Bgs. (29)-(3l) , 
are  physically  genuine  solutions  pertaining  to  the  particular  (m,  p,  T_  )  data  set 
imposed.  u 
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1.6  THE  NONEQUILIBRIUM  NITRAMINE  EVAPORATION  LAW 


To  facilitate  independent  calculation  of  the  burning  rate,  m, 
with  the  present  model,  an  additional  physical  constraint  is 
required.  This  physical  condition  is  provided  by  the  nonequilibrium 
evaporation  law  at  the  melt/gas  interface. 

The  process  of  evaporation,  at  any  instant,  involves  molecules 
of  the  substance  (R)  leaving  the  liquid  in  an  outward  flux,  mR(out), 
and  an  influx  of  molecules  effectively  returning  to  the  surface, 
n^(in).  When  the  net  effect  is  n^fout)  -  ir^fin),  the  vaporization 
is  termed  at  equilibrium.  However,  when  m^(out)  *  m^fin) ,  net 
evaporation  takes  place,  as  it  is  assumed  to  occur  at  the  surface 
of  deflagrating  nitramines.  In  this  instance,  the  net  vaporizing 
material  flux  can  be  expressed  [37]  as, 

^  =  **PsUms[y^(o+)-Vrco*)]  <«> 

where  0  <  1  denotes  the  Knudsen  accomodation  coefficient,  ^s 

is  the  gas  density  at  the  surface?  the  mean  molecular  velocity 
perpendicular  to  the  surface  is 

^s  =  C^T(/2Trwe)'A  (33) 

eq 

and  Y  and  YR  denote  the  equilibrium-  and  actual-mass  fractions 
of  *R^  on  the  gas  side  of  the  surface,  respectively.  The  outward 
flux,  m^tout)  is  always  given  by  the  equilibrium  expression, 


=  [feVP]  Wt/Ws 


Wg  is  the  mean  molecular  weight  of  gas  at  the  surface,  and  pR  is  the 
partial  pressure  of  ,RI ,  where  use  have  been  made  of  Dalton's  law. 
Using  a  Clansius-Clapeyron  expression, 


=Tv  exp  [  J  ] 


(35) 


with  and  T^  being  the  measured  low  temperature  reference  values. 

The  net  vaporizing  flux,  therefore,  is  the  difference  between 

two  large  numbers,  since  U  U  .  Thus,  Y^g  and  Y  (actual)  must 

ms  gs  R  R 
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tit  comparable  to  obtain  the  moderate  values  of  m  observed;  indeed, 
Y  (0  |  is  expected  to  yield  a  good  estimate  of  the  actual  value 
or  Yr  (0  )  ,  but  their  difference  becomes  important  whenever  m  is  to 
be  calculated. 

Note  that  the  total  mass  flux  in  the  present  system  is 

vn  ~  vn]j  -Ye(cr)l 

accounting  for  subsurface  decomposition  in  the  melt  layer;  hence, 

'tn  =  vrvAe.Co')  (36) 

which  serves,  along  with  the  foregoing  formulation  leading  to 
as  the  auxiliary  condition  sought  for  independent  mass  flux 
calculations. 
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1.7 


Discussion  of  Preliminary  Results 


At  the  time  of  this  writing,  only  independent  melt  phase  results  are 
available,  presented  in  Tables  2  and  3.  Comparison  can  be  made  between 
solutions  which  include  bubbles  (Table  2)  and  those  which  exclude  bubbles 
(Table  3,  where  gaseous  decomposition  products  are  assumed  dissolved  in 
the  liquid)  at  three  distinct  pressures.  This  comparison  clearly  shows 
that  for  those  cases  which  include  bubble  formation:  (l)  The  extent  of 
subsurface  decomposition  is  appreciably  smaller.  (2)  The  heat  feedback 
at  the  surface  is  appreciably  larger  -  within  the  same  range  of  surface 
temperatures.  Overall  melt  thickness  obtained  in  both  bases  was  comparable 
for  the  same  p  and  Ts-range.  A  parametric  study  was  run  at  different  values 
of  T0  as  well;  the  results  are  withheld  due  to  uncertainties  in  the  corre¬ 
sponding  values  of  the  burning  rates,  which  have  to  be  imposed  (using 
experimental  data). 

Several  numerical  methods  were  attempted  for  the  generation  of 
independent  gas  phase  solution  manifolds.  The  quasi-nonsteady  approach 
which  was  previously  used  with  success  [23-25]  (in  which  a  parabolic  partial 
differential  system  is  solved  by  marching  in  psuedo-tirae  toward  steady  state) 
was  found  too  cumbersome  in  terms  of  core  requirements  and  CPU  time  and 
tended  at  times  to  diverge  with  the  diverse  kinetics  data  employed.  Direct 
initial-value  methods  (marching  forth  in  space  from  the  liquid/gas  surface), 
which  are  much  less  costly,  are  intrinsically  divergent,  owing  to  the  presence 
of  one  set  of  eigensolutions  which  diverge  as  large  distance  from  the  sur¬ 
face  is  attained;  this  occurs  for  the  logarithmic-transformed  system,  Eq. 

(26) ,  as  well. 

At  present,  a  transform  to  the  4 ?  -1  plane  has  been  performed  to 
use  in  a  direct  (low-cost)  initial-value  solution  for  the  gas  phase  system 
in  an  iterative  mode.  This  approach  is  fully  formulated,  and  is  expected  to 
overcome  the  previous  numerical  difficulties  associated  with  the  diverse 
chemical  kinetics  scales  imposed. 


57 


References 


I. .  Robertson,  A.  J.  B.,  "The  Thermal  Decompositon  of  Explosives, 

Part  II.  Cyclotrimethylenetrinitramine  and  Cyclotetramethylenetetrimitramine 
and  Cyclotetramethylenetetranitramine,"  Trans.  Faraday  Society,  No.  1+5, 
pp.  85-93,  19^9 • 

2.  Rauch,  F.  C.  and  Fanelli,  A.  J.  ,  "The  Thermal  Decomposition  Kinetics  of 
Hexahydro,  1,  3,  5-trinitro-s-triazine  Above  the  Melting  Point: 

Evidence  for  Both  a  Gas  and  Liquid  Phase  Decomposition/'  The  Journal 
of  Physical  Chemistry,  Vol.  73,  No.  5,  May  1969,  pp.  160U-1608. 

3-  Rogers,  R.  N.  and  Smith,  L.  C.  ,  "Application  of  Scanning  Calorimetry 
to  the  Study  of  Chemical  Kinetics, "  Thermochimica  Acta,  No.  1,  1970, 
pp.  1-9. 

b.  Rogers,  R.N.,  "Differential  Scanning  Calorimetric  Determination  of 

Kinetic  Constants  of  Systems  that  Melt  with  Decomposition,"  Thermochimica 
Acta  3,  1972,  pp.  1437-^7.  - 

5.  Roger,  R.  N.  and  Daub,  G.  W. ,  "Scanning  Calorimetric  Determination 

of  Vapor-Phase  Kinetic  Data,"  Analytical  Chemistry,  Vol.  b 5  pp.  596-600, 

March  1973. 

6.  Rogers,  R.  N.,  "Determination  of  Condensed  Phase  Kinetic  Constants," 

Note,  Thermochimica  Acta,  No.  9,  PP,  UM-446,  197*+. 

7*  Suryanarayana,  B. ,  Graybush,  R.  J.  and  Autera,  J.  R.  "Thermal  Degradation 
of  Secondary  Nitramines :  A  Nitrogen-1  Tracer  Study  of  HMX  (l,  3,  5,  7 
Tetranitro-1 ,  3,  5,  7-  Tetraazacyclooctane^ "  Chemistry  and  Industry,  pp. 
2177-2178,  December  30,  1967. 

8.  McCarty,  K.  P. ,  "HMX  Propellant  Combustion  Studies:  Phase  1,  Literature 
Search  and  Data  Assessment  AFRPL-TR-76-59 ,  December,  1976. 

9*  Cosgrove,  J.  D.  and  Owen,  A.  J.,  "The  Thermal  Decomposition  of  1,  3,  5 
Trinitrohexadydro-1 ,  3,  5  triazine  (RDX), "  Chemical  Communications, 

1968,  p.  286. 

10.  Cosgrove,  J.  D.  and  Owen,  A.  J.,  "The  Thermal  Decomposition  of  1,  3, 

5  Trinitro  Hexahydro  1,  3,  5  Triazine  (RDX)  Part  I  The  Products  and 
Psysical  Parameters  Combustion  and  Flame,  No.  22,  197U,  pp.  13-18. 

II.  Cosgrove,  J.  D.  and  Owen,  A.  J. ,  "The  Thermal  Decomposition  of  1, 

3,  5  Trinitro  Hexahydro  1,  3,  5  Triazine  (RDX)  Part  II:  The  Effects 
of  the  Products  ,M  Combustion  and  Flame,  No.  22,  197I+,  pp.  19-22. 


Rogers,  R.  N.  and  Morris,  E.  D. ,  "On  Estimating  Activations  Energies 
with  a  Differential  Scanning  Calorimeter  , 11  Analitical  Chemistry,  Vol. 
38  (19 66),  pp.  !|12-Ul4. 


58 


13.  Pollard,  F.  H.  and  Wyatt,  R.  M.  H. ,  "Reactions  Between  Formaldehyde 
and  Nitrogen  Dioxide.  Part  I:  The  Kinetics  of  the  Slow  Reaction," 

Trans .  Faraday  Society,  Vol.  1*5,  191+9,  pp.  76O-767. 

ll+.  Pollard,  F.  H.  and  Woodward,  P.  ,  "Reactions  Between  Formaldehyde 
and  Nitrogen  Dioxide.  Part  II:  The  Explosive  Reaction,"  Trans . 

Faraday  Society  45,  19I+9 ,  pp.  767-770. 

15.  Pollard,  F.  H.  and  Wyatt,  R.  M.  H. ,  "Reactions  Between  Formaldehyde 
and  Nitrogen  Dioxide.  Part  III:  The  Determination  of  Flame  Speeds," 
Trans .  Faraday  Society  1+6,  1950,  pp.  281-289. 

16.  Goshgarian,  B.  B. ,  "The  Decompositon  of  Cyclotrimethylenetrinitramine 
(RDX)  and  Cyclotetramethylenetetranitramine, "  AFRPL-TR-78-76 ,  October 
1978. 

17.  Brenecker,  R.  R.  and  Smith,  L.  C. ,  "On  the  Products  Formed  in  the 
Combustion  of  Explosives.  Freeze  Out  of  the  Water-Gas  Reaction," 

Journal  of  Physical  Chemistry,  71,  No.  8,  1967,  pp.  2381-2390. 

18.  Rocchio,  J.  J.  and  Juhasz,  A.  A.,  "HMX  Thermal  Decomposition  Chemistry 
and  Its  Relation  to  HMX-Composite  Propellant  Combustion,"  11th  JANNAF 
Combustion  Meeting,  CPIA  Publication  No.  26l,  December  197I,  pp.  2kj- 
265. 

19.  Lenchitz,  C.  and  Velicky,  R.  W. ,  "The  Role  of  Thermochemistry  in  HMX 
Propellant  Burning,"  12th  JANNAF  Combustion  Meeting,  CPIA  Publication 
273,  Vol.  II,  Dec.  1975.  PP-  301-321. 

20.  Shaw,  R.  and  Walker,  F.  E. ,  "Estimated  Kinetics  and  Thermochemistry 
of  Some  Initial  Unimolecular  Reactions  in  the  Thermal  Decompositon 
of  1,  3,  5,  7  Tetranitro-1 ,  3,  5,  7  tetraazacyclooctane  in  the  Gas 
Phase,"  The  Journal  of  Physical  Chemistry,  Vol.  8l,  No.  25.  no.  2572- 
2576,  1977. 

21.  Fifer,  R.  ,  "Shock  Tube  Study  of  the  High  Temperature  Kinetics  and 
Mechanism  of  Nitrogen  Dioxide  -  Aldehyde  Reactions , "Modern  Development 

Shock  Tube  Research ,  Proc.  10th  International  Shock  tube  Symposium, 
Kyoto,  Japan,  1975,  pp.  613-620. 

22.  Fogelzang,  A.  E. ,  Svetlov,  B.  S. ,  Azhemian,  V.  Ya.  Kolyasov,  S. 

M.  and  Sergienko,  0.  I.,  "The  Combustion  of  Nitramines  and  Nitosamines, " 
Translated  from  Docladi  Academi  Nauk,  SSSR,  Bol.  216,  No.  3  May  1971, 
pp.  603-606. 

23.  BenReuven,  M. ,  "Nitramine  Monopropellant  Deflagration  and  Nonsteady, 
Reacting  Rocket  Chamber  Flows,"  Ph.D.  Thesis,  No.  ll+55-T  Dept  of 
Mechanical  and  Aerospace  Engineering,  Princeton  University,  Princeton, 

NJ  November  1979. 

2b.  BenReuven,  M. ,  Caveny.  L.H.,  Vichnevetsky ,  R,  and  Summerfield,  M.  , 

"Flame  Zone  and  Subsurface  Reaction  Model  for  Deflagrating  RDX," 
Proceedings,  l6th  Symposium  (International)  on  Combustion,  The  Combustion 
Institute,  Pittsburgh,  PA,  1976,  pp  1223-1233. 


59 


25-  BenReuven,  M.  and  Caveny,  L.H.,  "Nitramine  Flame  Chemistry  and 

Deflagration  Interpreted  in  Terms  of  a  Flame  Model, "AIAA  Jour.  Vol. 

19,  No.  10,  Oct.  1981,  pp.  1276-1285. 

26.  Schroeder,  M.A.  "Critical  Analysis  of  Nitramine  Decomposition  Results: 

Some  Comments  on  Chemical  Mechanism,  "Proc .  l6th  JANNAF  Combustion 
Meeting,  CPIA  Pub.  No.  308,  Vol.  II,  Dec.  1979,  pp.  17-31* 

27.  Schroeder,  M.A.  "Critical  Analysis  of  Nitramine  Decomposition  Data: 
Activation  Energies  and  Frequency  Factors  for  HMX  and  RDX,"Proc.  17th 
JANNAF  Combustion  Meeting,  CPIA  Pub.  329,  Sep.  1980,  Vol.  II,  pp.  1*03- 
508. 

28.  Schroeder,  M.A.  "Critical  Analysis  of  Nitramine  Decomposition  Data: 

Product  Distributions  from  HMX  and  RDX,"Proc.  JANNAF  Combustion  Meeting, 
CPIA  Pub.  3l*7,  Oct.  1981,  Vol.  II,  pp.  395-1*13. 

29.  McGuire,  R.R.,  and  Traver,  C.M.  "Chemical  Decomposition  Models  for 
the  Thermal  Explosion  of  Confined  HMX,  TATB,  RDX,  and  TNT  Explosives," 
Proc.  7th  Detonation  Symposium  (international),  Annalpolis  MD,  June 
1981,  pp.  550-557. 

30.  Beyer,  R.A. ,  and  Morgan,  C.U.  "Electron  Spin  Resonance  Studies  of 
HMX  and  RDX  Decomposition, "Proc .  l6th  JANNAF  Combustion  Meeting,  CPIA 
Pub.  308,  Dec.  1979,  Vol.  II,  pp.  51-57 

also,  cf 

Pough,  H.L.  Davis,  L.P. ,  Wilkes,  J.S.  ,  Carper,  W.R.,  and  Dorey,  R.C. 
"Thermal  Decomposition  of  RDX  Below  the  Melting  Point, "Proc.  7th 
Symposium  (international)  on  Detonation,  Annapolis,  MD,  June  1981, 

PP-  1*5-52. 
also,  cf 

Axworthy,  A.E.,  Flanagan,  J.E.,  Woolery,  D.O.  "High  Temperature 
Pyrolysis  Studies  of  HMX,  RDX  and  TAGN, "Proceedings ,  15th  JANNAF 
Combustion  Meeting,  CPIA  Pub.  297,  Sept.  1978,  pp.  253-265. 

31  Farber,  M.  and  Srivastava,  R.D.  "Thermal  Decomposition  of  HMX," 

Proceedings,  l6th  JANNAF  Combustion  Meeting,  CPIA  Publication  308 
Dec.  1979  Vol.  II,  pp.  59-71. 

32.  Fifter,  R.A.  and  Holmes,  H.E. ,  "Kinetics  of  Nitramine  Flame  Reactions," 
Proceedings,  l6th  JANNAF  Combustion  Meeting,  CPIA  Pub.  308  September, 

1979,  Vol.  II,  pp.  35-50. 

33.  Derr,  R.  L. ,  Boggs,  T.  L. ,  Zurn,  D.  E.  and  Dibble,  E.  J. ,  "The 
Combustion  Characteristics  of  HMX,"  11th  JANNAF  Combustion  Meeting, 

CPIA  Publication  261,  Vol  1,  Dec.  197h,  pp.  231-21*1. 

3l*.  Nigmatulin,  R.I.  "Methods  of  Mechanics  of  a  Continuous  Medium  for 

the  Description  of  Multiphase  Mixtures , "PMM  Vol.  3I*,  No.  6,  1970,  pp. 
1097-1112  (TRANS. ) 

35.  Williams,  F.  A.,  Combustion  Theory,  Addison-Wesley ,  Reading,  Mass., 

1965,  Chapter  1:  "Summary  of  Basic  Fluid  Dynamics  and  Chemical  Kinetics," 
PP-  1-17- 


60 


36.  Scala,  S.  M.  and  Sutton,  G.W. ,  "Energy  Transfer  at  a  Chemically 
Reacting  or  Slip  Interface,"  ARS  Journal,  Fet.  1959*  pp.  ll+l-1^3* 

37*  Plesset,  M.S.  and  Prosperetti ,  A. ,  "Bubble  Dynamics  and  Cavitation," 
Annual  Reviews,  Fluid  Mechanics,  9»  1977s  PP*  1^+5— G 8 5 - 


61 


TABLE  1 


o 

w 


oo 


o 

rH 

X 

co 


in 


ip 


' D 

m 

rH 

2 


* 

,<U 

MH 

•H 

pH 


♦ 

♦ 

U 


£ 


0J 


2 


in 

• 

in 

• 

• 

• 

* 

a\ 

• 

00 

• 

r- 

iH 

<Ti 

r- 

O 

CN 

00 

rH 

CN 

r** 

in 

in 

rH 

rH 

CO 

rH 

U 

O 

O 

rH 

O 

0  4J 

3 

in 

3 

• 

co 

o 

O 

3 

3 

00 

in 

• 

ip 

o 

• 

3 

in 

» 

IP 

• 

co 

rH 

00 

CN 

rH 

CN 

rH  C 

(N 


CM 


CN 


CN 


CN 


CN 


§ 

4* 

CN 

2 

4- 

§ 

4- 

o 

CN 

w 

CN 

O 

O 

O 

4- 

+ 

CN 

B 

ef 

B^ 

8 

CN 

8 

CN 

ro 

5* 

in 

Cp 

62 


TABLE  2. a  RDX  MELT  WITH  GAS  BUBBLES 
Tq  =  300  K,  p  =  2  MPA  m  =  8.76  kg/m2-s 
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9  CEh 

21 

0 . 66 73 E 4  03 

c. 

8GCCE4  0  0 

0. 9991E4  00 

0. 97C4E4 01 

0. 6453E4  07 

C .  1 6  IJ  4  E  h 

a  cl 

0 .  6769E+ 03 

0. 

84  0  0E4 0  0 

0. 99 84 E 4  0  0 

0 .101 5E4  C£ 

:  0. 5988E+ 07 

C .  1 4 

72Eh 

22 

0. 6864E+  03 

0. 

88  00E4  0  0 

0. 997 1E4  00 

O. 1 C64E+ 0c 

:  0. 5££7E+  07 

C .  1 2 

:77Eh 

24 

C.6958E4C3 

c. 

9ECCE4  00 

0 . 994 6E 4  o  0 

0. 1 1£3E4 0c 

:  0.411 1E4  07 

C .  1  C 

:C6Eh 

8w 

0.  7  C53E-+  03 

0. 

96CCE4  00 

0. 9886E4  00 

C .  1 2 C6Eh  02  C .  266 1  Eh-  07 

C  •  6  6 

-39  Eh 

26 

C.7147E+C3 

c. 

1 CCCE4  01 

0. 9648E4  0  0 

0 .  1 366 E+ 0£ 

:  0. 1 C£8E4 07 

U .  c'fci 

131  Eh 

KTDT  ?  J.-"M£-  S 

I:  TR  U.-"  I:  Z 

U6  j  M.--S 

PS  I  -  VG-- "V 

1  KE^KCLIQ) 

i 

0 .  8 C78E4 07 

c . 

1 1 17E4  01 

0. 5CC3E4 00 

c .  c 

0. 1 CC0E4  01 

c! 

0 . 8S34E+ 07 

c . 

1 157E4  01 

0.51 C£E+ 0  0 

c .  c 

0. 1 0  0CE4  01 

W 

0. 839 CE 4 07 

0. 

1197E4C1 

0. 5£01E4  00 

c .  c 

0. 1 CCCE4  01 

4 

C.854  6E4C7 

c. 

1  £  3  7  E  4  0 1 

0. 5£99E4 00 

c .  c 

0. 1 CCCE4  01 

cr 

0 . 8  7  0 1 E  4  0  7 

0. 

1E77E4C1 

0. 53 98 E 4 00 

c .  c 

0. 1 CCCE4 01 

6 

0 . 8 8  5  7 E  4  o  7 

c. 

1317E4  01 

0 . 54  97E4  0  0 

c .  c 

0. 1 C00E4  01 

r' 

C.9C13E4C7 

c. 

1357E4  01 

0. 5596E4 00 

c .  c 

0. 1 CCCE4 01 

8 

0. 9169E4  07 

c. 

1397E4C1 

0. 5695E4 00 

IJ  -  IJ 

0. 1 CCCE4 01 

9 

0.  9384 E4  07 

c. 

1437E+ 01 

0. 5793E4  0  0 

C  -  c 

0. 1  CODE 4 01 

1  C 

0 . 94  79E4 07 

0. 

14  76E4  01 

0 . 589£E4  0  0 

C. 1163E-C3 

:  0. 9 9 99 E 4  00 

1 1 

0 . 9634  E4  07 

0  - 

1 5 1 6E4  0 1 

0. 599 1E4 uC 

C. 2372E- C2 

:  0. 99 98 E 4 00 

12 

0. 9787E4 07 

c. 

1555E4  01 

0. 6C9CE4 00 

C. 6U23E- C3 

:  C.9995E+CC 

13 

0. 9938E4  07 

C. 

1594E4  01 

0. 61S8E4 00 

C.  1 1 C2E- 02 

:  0. 9991E+ 00 

14 

0. 1 0C8E4 08 

C. 

163EE4C1 

0 • 6£87E4 0  0 

C .  2  2  -  •  7  E  “  C2 

:  C.9981E4C0 

13 

0.  1 C££E4 08 

c. 

1667E4C1 

0. 6386E4  0  0 

0 « 4  284  E- 02 

:  0. 9964 E 4 00 

16 

C.  1C35E4C8 

0. 

170  UE4 C 1 

0. 6485E4 00 

C. 7778E- 02 

:  0. 9935E4  00 

17 

0. 1 C45E4  08 

C. 

1 7E6E+  0 1 

0. 6584E4  0  0 

C.  14  CEE- Cl 

0 . 9’  8  8  £  E  4  0  0 

18 

0.  1  C58E4  08 

c. 

1743E4C1 

0. 6682E4 00 

C .  ciA  67E-  C 1 

0. 9 7 93 E 4  00 

19 

0.  1 C58E4  08 

C. 

1744E4  01 

0 . 678 1 E4  0  0 

C.4259E-C1 

C.964EE4 00 

2C 

0.  1 C48E4  08 

0. 

1719E4  01 

0. 688 CE 4 00 

C. 721  IE- Cl 

0. 9394 E 4 00 

21 

0 .  1 0 1 8E4 08 

0. 

1656E4  01 

0. 6979E4  0  0 

C. 1 194 Eh CC 

0 . 8997E+ 00 

22 

0.  97C9E4  07 

0. 

1537E4  01 

0. 7078E4  00 

C. 1923E4CC 

0. 838£E+ 00 

£2 

0  ■  9 4  3E4  07 

C. 

134  8E4C1 

0 . 7 1 7  6  E  4  ij  0 

C. 3003 Eh CC 

0.  74  74E4  0  0 

24 

0 .  78 1 9E4 07 

c. 

1C55E4C1 

0. 7875E4  0  0 

C.45C2E+CC 

0 .  6£ 1 5E+  0 0 

23 

0 . 634 6E4 07 

0. 

683CE4  00 

0. 7374E4 00 

U -  6394 Eh  C C 

0. 46£5E4 00 

26 

0. 46£5E4 07 

0. 

£64  CE4 0  0 

0. 7473E4 00 

C.  84 99 Eh  CC 

C.£855E4CC 

63 


~4 

A 

A 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


TABLE  2.b  RDX  MELT  WITH  GAS  BUBBLES 


T0  =  300  K,  p  =  6  MPa,  m  =  19.7  kg/m2-s 


4 4 CO HE EH SEE  PURSE  PRDFILES+4 


1 

TEMP  ao 

DMLSS. TEMP 

VMTR 

V »  M I  CRD  NS  KET/'EY  >  J--M8- : 

S  RHDE KC--  M3 

1 

c 

. 4  785E4  03 

0  -  0 

0. 1 0CGE 4 01 

0-  0  . 

0. 9798E4 07 

0. 188CE+ 04 

lL 

c 

.  4  894  E 4  03 

C. 4  OCCE- 01 

0. 1 CCCE4 01 

0  *  3 1 83 E+  0  0 

0. 1 C8  0E4  08 

0. 188CE4 04 

c 

. 5CC8E4  03 

0.80 COE- 01 

C. 1 OCCE 4 01 

C.  6  £  43  E  C  0 

0.1 06  CE4 08 

0. 188CE4  04 

4 

c 

.  51 1 1E4 03 

0.  18 COE 4 00 

C. 1 OCCE 4  Cl 

C.  9189E+  CC 

0. 11 0CE4  08 

0 .  1 88 CE4 04 

cr 

L 

.  5819E4  03 

0.  16CCE4  00 

0. 1  CODE 4  ci 

0 .  I  £  C3E-+  0 1 

0.  114CE4C8 

0 .  1 88  CE4 04 

t 

c 

. 5388 E 4 03 

0. 3CCCE4 00 

0. 1 OCCE 4 01 

C  •  14  77E-+  0 1 

0.  1181E4C8 

0.  188CE+C4 

i*' 

f 

. 5437E4 03 

0.84  CCE4  00 

0. 1 CCCE4 01 

0 . 1742E+C1 

0 • 1 £8 1 E4  08 

0. 188CE4 04 

8 

r 

.  334  3E-+  C3 

C .  £8  0  CE+  0  0 

0.  1 OCCE 4 01 

G. 1998E-+ 01 

0 .  1 8 6 1 E4  08 

0.  188CE4  04 

o 

c 

. 5654E4  03 

0. 38 COE 4 00 

0 .  1  0  0  CE4  0 1 

C . ££46E+  C 1 

0.  1 3 C8E4  08 

0 .  1 88  OE4  04 

1  C 

c 

. 576EE+  C3 

0. 36CGE+ 00 

0. 1 OCCE 4 01 

0 .  £ 4 8 7 E  +  C 1 

0.  1348E4  08 

0.  188CE4 04 

1 1 

L 

.  587 1 E4  03 

0. 4CCCE4  00 

0. 1 OCCE 4 01 

C.  2781E-+  01 

0 . 1 388E4 08 

0.  18ECE+ 04 

1£ 

f 

. 598 CE4 03 

0. 44CCE4 00 

0. 1 OCCE 4  01 

C.  £947E*+  Cl 

0 . 1 488E4 08 

0. 1819E4  04 

IS 

c 

. 6C88E4  03 

0. 48CCE4 00 

0.  1 CCCE4  01 

0 . 3 1 68E+ 0 1 

0 . 1 468E4 08 

0. 1819E4  04 

14 

f 

-  6 1 87E+ 03 

0. 52CGE-+  00 

0.  1  CODE 4  01 

G . 3 3 83 E 4 G 1 

0. 15CCE4 08 

0. 1817E404 

15 

c 

. 63C5E4  03 

C. 56CCE+ 00 

C • 9999E+ G G 

C.3592E+C1 

0.  1538E4C8 

0. 1815E4  04 

16 

c 

.6414E4C3 

0. 6 CODE* 00 

0. 9999E4  00 

0  •  3797E-+  U 1 

0.  1573E4C8 

0.  181CE4C4 

17 

c 

. 63£3E+ 03 

0. 6-4CCE4  00 

C - 9998E+ G G 

C .  3997E-+  C 1 

0.  16C4E4C8 

0 .  1 8  0 1 E4  04 

18 

f 

.  663 IE-*-  03 

C.  68 COE**  00 

0 . 9996E4  0 0 

C.4194E4C1 

0 . 1 688E4 08 

0 .  1 786E4 04 

18 

r 

.  6  7  4  U  E  ■+  IJ  3 

0. 78CCE4 00 

0. 9 9 93 E 4  00 

0 . 4 3 8 8 E+  Cl 

0. 1639E4 08 

0. 1759E4 04 

£  C 

c 

- 6848E+ 03 

0. 76CGE+ 00 

C. 9987E4  00 

C.4S83E4C1 

0. 163CE4  08 

0. 1718E4  04 

£1 

. 6957E4  03 

C. SOCCER  CO 

C - 9977E+ C C 

0 •  47 8 C E-*-  Cl 

0.  1588E4C8 

0. 1636E4  04 

lL  l! 

.  7  066E4 03 

0. 84  COE 4 0  0 

C .  996  GE+  G  C 

C. 4985E+C1 

0.  1498E4C8 

0. 1514E4  04 

C. 

. 7174E4  03 

0. 88CCE4 00 

Cm 9988E+ GO 

0.S2C9E4C1 

0 . 1 348E+ 08 

0. 1331E4  04 

£4 

. 7883E4 03 

0 . 98  0 CE4  0  0 

C.  9865E-*  0  0 

C. 5468E4 01 

0.  11 C6E4  08 

0. 1 C78E4  04 

•”  cr 

I _ W 

i 

.7391E4C3 

0.  96CCE4  0  0 

0 • 97£6E+ 0 C 

C- S8C3E4 01 

0 . 798 1 E4  07 

0 . 7 4  6 7 E  +  03 

£6 

i 

. 75 COE 4  03 

0 .  1  C  0  0E4  c  1 

0  -  9333E-+  00 

0 . 63 1 9  E *+•  0 1 

0. 46S8E4 07 

0.4  0 04 E 4 03 

i 

HTCTj J/K2-S 

DTRIVEZ 

U6 N/S 

psi-ve^v 

KE.-  K  <L  I  C> 

1 

£ 

1818E4  08 

0 . 9 7 1 6E4 0 C 

“oT^TOTToiT 

Cm  C 

0.  1  0 cOt-4  Ci 

c. 

( 

1859E4  08 

0. 1 018E4  01 

0. 3839E4 00 

0 .  0 

0. 1 OCCE 4  01 

( 

1899E4C8 

0.  1C58E4C1 

0 . 3984  E+  0 0 

0 .  0 

0. 1 OCCE 4  01 

4 

( 

J.  1939E4  08 

0.  1 098 E4  01 

0.4  CC9E4  00 

0.  0 

0. 1 OCCE 4 01 

cr 

i; 

1979E4  08 

0 .  11 3 8  E 4 0 1 

0. 4  094 E 4 0  0 

0 .  0 

0. 1 OCCE 4 01 

6 

( 

8C8CE4  08 

0.  1178E4C1 

0. 4 18CE4 00 

0 .  0 

0. 1 OCCE 4  01 

( 

8C6CE4  08 

0. 1818E4  01 

0. 4  865E4  00 

c.  0 

0. 1 OCCE 4 01 

e 

( 

5.81  0  0E4  08 

0.  1858E4  01 

0.4  35CE4C0 

0.  0 

0. 1 OCCE 4 01 

8 

( 

8141E4  08 

0. 1891E4  01 

0. 4  4  35 E 4  0  0 

0.  0 

0. 1 OCCE 4 01 

1C 

( 

5. 8181E4  08 

0.  1331E4C1 

0. 4  580E4  0  0 

C.7913E-C4 

0. 9999E4  Ci 

1  1 

( 

J »  c'c  c!  1  E  +  IJ  8 

0.  1371E4C1 

0. 46C5E4  00 

C. 1619E- 03 

0. 9999E+  Ci 

1£ 

i; 

J . ££6 1 E+ 08 

0 .  1 4  1 1 E4  0 1 

0. 4691E4  00 

0- 3300E- 03 

0 . 9997 E 4 0 i 

18 

( 

j.  £3  0 1 E  +  08 

0.  145CE4C1 

0 . 4776E+ 0 0 

|J  •  7 1 36E-  03 

0.  9994 E 4  Ci 

14 

i 

J . £339E+ 08 

0.  14  89E4C1 

0 . 4  86 1 E4  0  0 

C.  1493E-C2 

0 . 9 9 8 7 E  4  oi 

13 

i 

J .  2377E+  08 

0.  1586E4C1 

0. 4  94  6E4 00 

C .  £9 09 E- 0£ 

0. 99 76 E 4  oi 

lb 

i 

J .  £4  1  cE  +  08 

0. 1561E4  01 

0.5C31E40C 

0 . 5394 E- 0£ 

0.  9953E4  Ci 

17 

i 

;.£4  4SE4C8 

0.  1598E4C1 

0.51  17E4CC 

0.  104  8E-01 

0 ■ 99 1 8E4 0 i 

18 

i 

J . £4  66E+ 08 

0. 1615E4C1 

0 . 58  C8E4  0 0 

0. 1923E- 01 

0. 9838E4 Ci 

18 

i 

.!,  iz!  4  7  7  E  *+  U  fc! 

0 . 1 686E4 0 1 

0. 5887E4 00 

C.3447E-C1 

0. 971 CE4 Oi 

£C 

i 

I' .  £4  68E*+  08 

0 . 1 6 1 7E4  0 1 

0 . 5378E4 0  0 

0 . 6  039E- 0 1 

0. 94  9SE4  Ci 

£1 

i 

J.  £4£3E-fr  C8 

0.  1576E4C1 

0.54  57E4CC 

0 . 1 C2£E+ 0 C 

0. 9138E+  Ci 

££ 

i 

C .  £334  E-+  08 

0.  14  87E4  01 

0. 5543E4  00 

0. 1712E+  0  0 

0. 856 CE 4 0 

C.  w 

i 

C . £  1 73E4  08 

0. 1 338E4 01 

0. 5688E4  00 

0 .  £74  1 E-+  L*  C 

0. 7696E4 0 

£4 

i 

C.  1 934  E+  08 

0.  1 098 E 4  01 

0. 5713E+  0  0 

0.41 87E^ 0 U 

0 . 64  8 1 E4  0 

•“  cr 

i_ 

i 

0.  1614  E+  08 

0 . 79 1 9E+ 0 0 

0. 5 7 98 E 4 00 

C. 6CC9E+ 00 

0 ■ 4 948E4 0 

£b 

i 

0. 1£31E+  08 

U .  4  6  4  6  E  4  i_!  i_! 

0. 5 8 83 E 4  0  0 

0 . 7 9 4  6  E  4  C  IJ 

0. 338 OE 4 0 

64 


TABLE  2 ♦ c  RDX  MELT  WITH  GAS  BUBBLES 
Tq  =  300  K,  p  =  10  MPa,  m  =  30.2  kg/m2-s 


■mCDMBEMSEE  phrse  profiles-*--* 


1 

TEMP  00 

EMLSS. TEMP 

YMTR 

Y »  MI CRCMS  KET.-  EY «  J/K2- S 

RHCE  -  KG'-" M3 

1 

C.4  785E4C3 

c.  c 

c 

1 CCCE 4 01 

e.  c 

0. 15CCE4 08 

0. 18£CE4 04 

C 

C.4894E4C3 

0.4 CCCE- 01 

c 

1 CCCE4  01 

C. 2C78E4 CC 

0.  1562E4C8 

0. 188CE4  04 

p 

C. 5CC£E+ C3 

0. 80 COE- 01 

c 

1 CCCE 4 01 

C  •  4  C  7  6  E  4  C  C 

0. 1623E4 08 

0 .  1 8£  CE4  04 

4 

C.51UE+C3 

0. 12CCE4 00 

c 

1CCCE4C1 

C.5999E+CC 

0.  1 685E+ 08 

0.1 8£ 0E+ 04 

cr 

C.5£19EtC3 

0 . 16C0E+ GO 

c 

1 0C0E4  Cl 

C.7352E4CC 

0. 1747E4  08 

C. 182CE4 C4 

6 

C. 53£8E4  03 

0. 2GC0E+ 00 

c 

1 CCCE 4 01 

C . 964 dE 4 C C 

0.  18C9E4 08 

0. 182CE4 C4 

»'* 

C.  54  37E-*  03 

0.34  C  *_!  E*+  CC 

e 

1 CCCE 4 01 

C. 1137E4C1 

C. 1871E4C8 

C .  1 8  £  0  E  4  ij  4 

g 

C. 5545E+  C3 

0. 880CE-*  00 

c 

. 1 CCCE 4  Cl 

C. 13C4E4  Cl 

0 . 1 932E4 08 

0.  1 8£ CE4 C4 

o 

C.  5654E-+  03 

0.32  00E-+  00 

c 

.  1  C C CE4  C 1 

C .  1 4 6 6 E  4  i_:  1 

0.  1994E4C8 

0. 18£CE4 04 

1  c 

C.  576EE-*  03 

0 . 36  0  CE-*  0  0 

c 

1 CCCE4  01 

C . 1 634 E 4 C 1 

0. 2C56E4 08 

0. 18£CE4 04 

1 1 

C.5871E+C3 

0.  4CCCE-+  00 

c 

1 CCCE 4  01 

C . 1 776E4  C 1 

0.211 7E4 08 

0. 182CE4 04 

12 

C. 598 CE* C3 

0. 44CCE+ 00 

c 

1 CCCE4  01 

C. 1934E4C1 

0 .  £  1 79E4 08 

0. 18£CE4 04 

13 

0 .  6 C88E+ C3 

0. 4SCCE+ 00 

c 

1 CCCE4 01 

0 . d C68E+ C 1 

0.  ££4  0E+  08 

0. 1819E4  04 

14 

C .  6197E-* C3 

0. 58CCE-+  00 

c 

. 1 CCCE 4 01 

C.22C8E+ Cl 

C.23C1E4 08 

0.  1819E4C4 

15 

C.  63C5E-*  03 

0.  56CCE-*  00 

c 

. 1CCCE4C1 

C.2343E4C1 

0 . £36 1 E+  08 

0. 1818E4  04 

It 

0. 64  14E-*  03 

0. 6CCCE-*  0  0 

c 

. 9999E4  0  0 

C .  2 4 7 8 E  4  p  1 

0 . £4 1 9E+  08 

0. 1816E4C4 

17 

C.6523E+C3 

0. 64C0E-*  00 

c 

. 9998E4  00 

C -  26 C8E4 C 1 

0. 2475E4 08 

0 .101 £E4 u  4 

18 

0. 663 IE + 03 

0. 68CCE-*  00 

e 

. 9997E+ C C 

C.2733E4C1 

C.2526E4C8 

0.  18C6E4  04 

19 

C > 674 C E  4 C 3 

C.72CCE+CC 

c 

. 9995E4CC 

Cm  226 CE4 C 1 

C.2568E4 08 

0.  1794E4  04 

2  0 

0 . 684  8E+  03 

0.  76 0 CE+  0 0 

c 

. 9991E4 00 

C • 2983E+ C 1 

0. £597E4  08 

0 .  1 77 4  E4  04 

£1 

C.6957E4  03 

0. 8CCCE-*  00 

c 

. 9984E4 CC 

C.31C3E4C1 

0 . £6  C3E4 08 

0.  174 1E+  04 

d  d 

0.7  C66E+ 03 

0. 84  COE -*  00 

0 

. 9973E+ 0 0 

C. 3228E4 Cl 

C.2574E4 08 

0.  1686E4  04 

C.  w 

0 .  717 4  E-*  03 

0. 88 COE 4 00 

c 

. 9954 E4  00 

C. 3353E4 Cl 

0 .  £4  92 E 4  08 

0. 16 COE 4 04 

£4 

0. 7£S3E+ 03 

0. 92CCE4 00 

p 

. 9922E4 CC 

C.  34  83 E 4  C.1 

0 . 2 3 35 E 4 08 

0.  14  68E4  04 

£5 

0.  739 1E-*  03 

0 .  96  0  CE4  0 0 

e 

. 9863E4 0 0 

C. 3628E4 Cl 

0.  2 CSCE 4  08 

0 . 1 £77E4 04 

£6 

0.75  0  CE4  03 

0.  1  CCCE4 01 

c 

. 9752E4 00 

C  •  3 7 9 6  E  4  C 1 

0 . 171 7E4 08 

0.  1 C£1E4  04 

1 

KTDT * J/'ME-S 

ETRU-'EZ 

IJ6>  M--  S 

PS  I  ~  V6-*' V 

KE  -  K  <L I C> 

l 

C. 2785E4 08 

0 . 9716E4 00 

c 

.34^6^400 

0 .  0 

0. 1 0CCE4 oi 

c. 

C .  £847E-*  08 

0.  1C18E4  01 

c 

.3528E4CC 

0.  0 

0. 1 CCCE 4 01 

W 

0 .  £9  0 9 E  ■*  08 

0. 1 C5£E4  01 

L 

36'  0  6  E  +  0  0 

0 .  0 

0. 1 CCCE 4  01 

4 

0.  £97  CE-*  08 

0. 1 C9£E4  01 

c 

.3685E4CC 

0.  0 

0. 1 CCCE4 01 

cr 

0. 3C3£E+ 08 

0. 1132E4 01 

c 

3763E4  0 0 

0.  0 

0.  1  CCCE 4  Cl 

6 

0. 3C94E-*  08 

0. 1 17£E4  01 

c 

384 1E4  00 

0.  0 

0.  1  CCCE 4 01 

0. 3156E-+  08 

0. 1212E4 01 

c 

392 CE 4 00 

0.  0 

0.  1 CCCE+ 01 

0 

0 . 3£  1 7E-*  08 

0 . 1 £5£E4  0 1 

c 

39 98 E 4 0  0 

0.  0 

0.  1 CCCE 4  01 

o 

0. 3£79E+ 08 

0.  1292E4  CM 

r 

4  C76E4 0 C 

0.  0 

0. 1 CCCE 4 01 

1  c 

0 . 334  1 E-*  08 

0. 1331E4  01 

c 

. 4154E4  0  0 

C.2348E-C4 

0. 1 CCCE 4 01 

1 1 

0 . 34  C£E+ 08 

0.  1371E401 

c 

. 4233E4 00 

C.729CE- 04 

C.9999E4 CC 

1  £ 

0.  3464E-*  08 

0. 141 1E4 01 

L 

.4311E4C0 

0. 148 IE- 03 

0. 9999 E 4  CC 

13 

0. 3525E4  08 

0. 1451E4  01 

c 

. 4389E4 00 

0 . 3 027E- 03 

C.9997E4 CC 

14 

0.  3586E-*  08 

0. 149CE4  01 

c 

. 4468E4  00 

0 .  6 1 59 E- 03 

0. 9 9 95 E 4 00 

15 

0 . 3646E-*  08 

0. 15£9E4  01 

c 

.4546E4CC 

0 . 1 ££6E- 02 

C. 999 CE+ 00 

16 

0.  37  04E-*  08 

0.  1567E4  01 

c 

. 4624E4  CC 

0. 2386E- 02 

0. 998 OE 4 00 

17 

0  •  376  0E-*  08 

0.  16G3E4  01 

c 

. 47C3E4 00 

0. 4467E- 02 

0.  9 9 6 £  E  4 cC 

18 

0.381 1E+ 08 

0. 1636E+  01 

c 

. 478 1 E4  C C 

0 . 8£3 CE- 0£ 

0 . 993 1 E4  0  0 

19 

C.3853E+ 08 

0. 1663E+ 01 

c 

. 4859E+ 00 

0. 148CE- 01 

0. 9876E4 00 

£  C 

0. 3881E-+  08 

0 . 1 68£E+ 0 1 

r 

. 4937E4 00 

0 ■ £6  G6E- 0 1 

C.9781E4 00 

£1 

0. 3886 E 4  08 

0 . 1 686E4 0 1 

. 5  0 1 6E4 c 0 

0. 4501 E-  01 

0. 96££E4 00 

££ 

0. 3856E+ 08 

0.1 667E+ 0 1 

. 5C94E4  00 

C.7597E-C1 

0. 9361E4 00 

£3 

0 . 377£E-*  08 

0. 1614E4 01 

.5172E4CC 

0. 1251E4  00 

0. 8949E+  00 

£4 

0 . 36 1  CE-*  08 

0. 1512E4 01 

. 5£51E+  0  0 

0.20 0 CE4  0 0 

0 ■ 83 1 9E4 0 0 

£5 

0. 3347E4 08 

0.  1347E4  01 

.5329E4CC 

C.3C82E4CC 

0. 74  C9E+CC 

£6 

0 . £97  CE+ 08 

0.  U12E4  01 

i 

.  54  C7E4 00 

C.4531E4CC 

0 . 6 1 9 1 E4  0  0 
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TABLE  3. a  RDX  MELT  W/OUT  BUBBLES 
Tq  =  300  K,  p  =  2  Mpa ,  M  =  8,76  kg/xn2-s 


■+  -+  L  cr  :r.-Er  ; -2-EI  PHPSE  PROFILED 

TEMP  fK'i  DMLSS. TEMP  YMTP 


1 

C.  4 785E.4  C3 

c.  c 

C.  1  CCCE 4  i.  i 

c. 

C  ■  4  8  7  4  E  4  0  3 

C. 4 CCCE- Cl 

C.  1  CCCE 4  Cl. 

0.  4  3*.  2E+  03 

0. SCCCE- Cl 

C. 1 CCCE 4 C j 

cr 

C.3053E+C3 

0. 314 1E+ 03 

0. 120CE4 00 
0. 16CCE4 00 

C. 1 CCCE 4  Cl 
C. 1 CCCE 4  Cl 

r. 

0 . 5  c  1  CE  4  03 

0. 2C00E+ 00 

C.  1 CCCE-*-  Cl 

C.  57  1  9E+  03 

0.  £4  CCE-+  CC 

C.  1 CCCE4 Cl 

8 

0.  3<  08E-*  03 

C. E8CCE+ CC 

C.  1  CCCE 4  i.:  i 

c* 

0 .  54  872*+  0  3 

C. 38 COE 4 CC 

C. 9999E* CC 

1  c 

0.  5586E-+  03 

f.,  ?.£. f:f:E-t  CC 

C.  9999E-*  C 

i  3 

*  c. 

<*.  54  7*E4  03 

0. 5765E4 03 

0. 4  CCCE 4 CC 
0.44 CCE4 00 

C. 9998E4 CC 
0 . 99 96E4 0 V 

J.  1 

.  5  75^E4  03 

0 . 4800E4 00 

0  •  9 9 9 2 E  t  C  o 

14 

15 

C.  3.-4  3E+  03 

0. 6  032E4  03 

0. 52CCE4 00 
0. 56C0E4 CC 

C. 9986E4 0 0 
0. 9974 E 4  00 

i  6 
17 

o  ■  t 1 2 1 E  *+  03 

0 .  62  1  CE4  03 

C. 6CCCE4 CC 
C.  64  CCE4 CC 

C. 9953E4 00 
0 .  99 1 7E4  o.O 

IS 

0  .  6-  2  2  *5'  E  4  0  3 

C,  68CCE-+  CC 

0.9857E400 

1 9 

C . 6388E4 03 

C.  ?£CCE4  CC 

0 . 9755 E 4  00 

c_  0 

0  .64  7  {'  E  4  0 3 

C .  76  C  CE+  C  C 

0 . 9 587 E 4 o  0 

21 

0 . 6566E4  03 

C .  8 C  C  CE-+  C  C 

0. 9312E4  0  0 

c’  2 

0. 6655E+ 03 

C.  84CCE+  CC 

0 • 8869 E 4 o  0 

£3 

0. 674 4 E4  03 

C.  88CCE-*  CC 

0 . 8 1 62E4  0  0 

£4 

0*  6833E+ 03 

C.  S'£CC.E+CC 

C*  7  04  1 E4  C  0 

O  cr 

Ui  w 

2* 

0.6 8  2 2 E  -+X*  3  - * 

0.7011 E4 03 

C.  1CCCE4C1 

i^C.'^£;69E4c.}t 
5':Cv.£372E+  CC 

,J 

HTCT  !•  J-^ME-S 

DTRLVCZ 

1 

C  •  8  0  7  8  E  *+  0  7 

C. 1 185E4 Cl 

2 

0  ■  *i‘ 2  2 5 E  4  1.!  i' 

C.  1££5E4  Cl 

c* 

0 • 8372E4 07 

C. 1£63E+  Cl 

4 

0. 85 1 9E4 07 

C.  13CSE4C1 

cr 

0 .  6- 6- 6E4  C 7 

0. 1343E+ Cl 

6 

0 . 88 1 2E4  07 

C. 1383E4  Cl 

»*' 

C.8959E4C7 

C.  14E5E+C1 

8 

0 . 9 1  0  6  E  4  i.:  7 

C. 14  63E4  Cl 

9' 

0. 9 2 52 E 4 07 

C.  13C5E+  Cl 

1  0 

0 . 8395 E 4 07 

C. 1343E4  Cl 

1  1 

0.9544E4C7 

C. 1384E4 Cl 

12 

0 . 5689E4 07 

C.  1 6£4  E4  C 1 

IS 

0 . 983 1 E4  07 

C .  1 663 E4  C 1 

14 

0  •  5  9  7  0 E  4  0  7 

C.  17C1E4C1 

15 

0 • 101 0  E  4  0  8 

C.  1739E4C1 

16 

0 .  1  U  2  2  E  4 08 

C.  1774E4  Cl 

17 

0.  1 033 E 4  08 

C .18 C6E4 C 1 

18 

0.  1  C4CE4  08 

C. 183EE4 Cl 

19 

0. 104  2E4  08 

C.  1848E4C1 

2  0 

0. 1 C36E4 08 

C.  184  8E4  Cl 

21 

0. 1 018E4 08 

C . 1 8£6E4 C 1 

22 

0. 9781E4 07 

C.  1763E4C1 

23 

0 . 9  063 E 4 07 

C.  164CE4  Cl 

24 

0 • 784 1 E  4  0  7 

C. 14£1E4  Cl 

£3 

C.5822E407 

C.  1C31E4C1 

£6 

0 .  24  25E4  l.:7 

C.  4SC8E4  CC 

MI  CP CMS  K 

DT.--DV.JVME-. 

*. 

C.4  331E4  iV;- 

58  c'kE-u  f:C 

C.  *«97E4  C7 

1 :  6 CE4 Cl 

C.  ^  r  ^  J 1 4  C7 

I  7  i  3  E  4  01 

L‘ .  4  7  5 1  E  4  0  7 

2 2 4 9 E 4  Cl 

0 . 4  c'E  4  07 

c  77  Ob*4  o  1 

C. 5  C8  5E4 0  7 

2276E4 0 1 

0. 5232E4 07 

77 6  7 Cl 

.  o. 5378E 4 07 

4  2*  4  6  E  o  1 

0. 5525E407 

«  1  c!  b  ■+  01 

0. 5672E-+  07 

5  I  66E-*  Cl 

0. 5818E4 07 

56  09 E4  0  1 

C . 594 3E 4 07 

6  4  1  t  4  r  1 

0. 61 06 E 4 07 

6 4 6 3 E 4  Cl 

0 . 624  8E4  07 

6.  6*  7  6-  E  4  f..'  l 

0. 6385E4 07 

7  2  8 1 E  4  01 

0. 6514E4  07 

7 6  7 7 E  4  01 

0 . 663  CE4 07 

.8 C68E4 01 

C • 6726 E 4  07 

84  54  E4  01 

0 .  67  86 E  4  07 

8 8 38b 4 01 

0. 679 CE4  07 

9224 E  4 01 

0 . 67  04 E 4  0 7 

9619E4 01 

C .  6  4  7  5  E  4  0  7 

1 0C3E4 02 

0  •  6  0 2 1 E  4  07 

1.C5CE+  CS 

Cm  52 1 6  E4  07 

1 1  C7£+T£. 

.0.  3858E+07 

i 2  C6E4 02 

C. 1545E407 

66 


ro  ro  ro  r»_*  ro  ro  ro 


TABLE  3.b  RDX  MELT  W/OUT  BUBBLES 
Tq  =  300  K,  p  =  6  Mpa,  m  =  1.97  kg/m^-s 


+  4CDKEEKSED  PHRSE  PRCFILES44 


t 

TEMP  OO 

DMLSS. TEMP 

YMTR 

V » M I OR  CHS 

1 

C.4  78SE4C3 

0.  0 

0. 1 CC0E4 Cl 

0.  0 

0. 9792E4 0 

l’ 

ill. 

0 . 4 89 1 E4  03 

0 . 4  0  0  CE-  0 1 

C. 1 COCE4 01 

0.3181 E+ 0 0 

0. 1 C19E4 0 

5 

0. 4998E+  03 

0.  8  0  0  CE-  01 

0. 1 CCCE4 01 

0 . 6 1 22E4 0 0 

0. 1 C38E4 0 

3 

4 

0 . 3 1  04  E4  03 

0. 1SCCE4 00 

C. 1 CCCE4 01 

0. 901 3E 4 00 

0. 1 C98E4 0 

f 

er 

0. 3c 1 1E4 03 

0. 1 6CCE4  00 

0.  1  0  0  CE 4  01 

0. 1 18CE4 01 

0 .  1  1  *'  7  E  4  u 

r* 

6 

C.5317E-+  03 

0. 2CCCE4  00 

0. 1 CCCE+ 01 

0.  143CE4  01 

0.  1  177E4  0 

s 

r 

0. 3484E4  03 

0. 24 00E4  0  0 

0. 1 CCCE4 01 

0. 171 CE4 Cl 

0 .  1 8 1 6  E  4  0 

: 

£ 

0. 553 CE4 03 

0. £8 COE 4 00 

0. 1 CCCE4 01 

0 .  1 963E4  ij  1 

0 .  1  8  3  6  E  4  f 

3 

o 

0. 3 6 SEE + 03 

0. 32CCE4 00 

0. 1 CCCE+ 01 

0. 2207E4 01 

0.  1293E4  0 

5' 

I  0 

0. 374 3E4 03 

0. 3E  00E4  00 

0. 9999E4  00 

0. 2444E4 01 

0.  1 333 E 4  0 

i- 

1  1 

0. 384 9E 4  03 

0. 4  CC0E4 00 

0. 9999E4  00 

0. 2674 E 4  01 

0. 1374 £4 0 

z 

2  z. 

0. 39 36 E +  03 

0.4  4  0  0E4  00 

0. 999 7E4  00 

0. 8.8 98 E  +  01 

0. 14  1384  0 

S 

X 

0. 6  C62E4  03 

0 . 4  8  0  0E4  0  0 

0. 9994 E 4 00 

0 . 2 1 1 6  E  4  0 1 

0.  14  32E4  0 

z 

A  1-1 

0 .  6  I  6  9  E  ■+  0  2 

0. 52 COE 4 00 

0. 9937E4  00 

0. 3328E4 01 

0.  1  4  90E+  0 

- 

*!  =r 

0. 6£75E+ 03 

0.  36  0CE4  00 

0 . 9976E4  0 0 

0. 3334 E4 01 

0.  1328E4  0 

z'- 

1 1£ 

C .  62S 1 E -►  03 

0  ■  E-  0  0  U  E  4  o  iJ 

0. 9954 E 4  00 

0. 3736 E 4 01 

0. 1363E4 0 

Z- 

X  1 

ij  •  t  A  c’  c'  E  +  0  3 

0. 64  COE 4 00 

0 . 99 13E4 0  0 

0. 3933 E 4  01 

0.  1393E4 0 

z'- 

i  £ 

0. *594E+C3 

0. 6800E4  00 

0. 984 £E 4 00 

L*  ■  4  1  cl  r"  E  4  IJ  1 

0 .  1 6  8 1  E  4  Ij 

5 

1 9 

C. 67C1E+ 02 

0. 72CCE+  0  0 

0. 9716E4 00 

0. 43 1SE+ 01 

0. 1636E4  0 

? 

S'  i’- 

0.  6SC7E-+  03 

0. 760  0E4  00 

0. 93 COE 4 00 

0 . 4  5  09 E 4  0 1 

0.  1634E4C 

=■ 

ii.  i 

0. 6914E+ 03 

0. 8CCCE+ 00 

0. 913264 00 

0 . 4  7  0 1 E  4  ij  i 

0.  16C4E4  0 

Z  z 

0.  7  OS  CE+  03 

0. 84  COE 4 0  0 

0.  3 5 SI  E 4  0 0 

U . 4 2 93 Et 0 1 

0 .  1 3  8  6  E  4  ij 

■: 

Z 

0.71 26E+ 02 

C.88CCE4  00 

0. 75 1 7E4 00 

0-51 14E4 01 

0.  1374  £4  Oi 

£: 

24 

0*  72 22 E 4  03 

0 . 92  0  CE4 0  0 

0. 5S95E4  0  0 

0 . 5364 E4  0 1 

0.11  03 E 4  Ci 

8 

•Z*  cr 

0 . 7339E+ 03 

0. 96C0E4  00 

0  •  3  3  4  cl  E  4  i ij  ij 

0. 5719E4  01 

0. 634 8E4 0 

r* 

S6 

0. 7446E+ 03 

0. 1 0CCE4  01 

0. 429CE- 01 

0.74 02 E 4 0 1 

0. 1373E4 0 

i' 

4 


l  u 


l  C 

X 

14 

15 

16 
1  7 
IS 
1 9 


HTDT*  J>T,£-S 

“  lSISE+CS 

C. 1658E4C8 
C .  1 8 9 7 E  4  C  8 
C . 1 937E+ C£ 
0 .  1  2  7  6  E  4  C  t' 
0 .  £  0 1 6  E  +  C  8 
C.2C5SE+C8 
C.  2093E-+  OS 
0 .  cl  z‘4E  +  OS 
0  •  2 1  7  4  E  +  0  8 
0 .SSI 2E+ OS 
0 • c  c! 5c'E4 0 c 
IJ «  2  c!  1  E  4  U 

0 .  S'  3  S  8  E  4  0  8 
0 .  S’ 3 6  5  E  4  OS 
0 > 2398E+ CS 
0  •  24  27E4  08 
0.  £*4  46E+  08 
0.  2451E4  08 
0. S421E+C8 
0 i 237CE+ 08 
0. £24 1E  + CS 
0.  SCME+  08 
0. 1598E+C8 
C.925SE4C7 
0. 173EE4C7 


DTRLVTZ 
C .  99  HE-t  0  0 
C .  1  03 1 E4  0 1 
C.  1C71E+C1 
0.11  UE4C1 
0.  U51E4C1 
0. 1 191E4 01 
0 . 1331 E+ C 1 
0. 1S71E4  Cl 
0.  131  1E4C1 
0. 1351E4 01 
0. 1391E4 01 
0. 1431E4 01 
0.  147CE4  01 
0.  13C9E+  01 
0. 1347E4  01 
0.  1582E4C1 
0.  1615E4C1 
0.  1641E4C1 
C.  1656E4C1 
0. 1633E4 01 
0.1  684  E+  0 1 
0.  1545E4C1 
0. 1391E+ 01 
0.  1117E+01 
0  ■  6-  6-  £'  9  E  4  ij  ij 
0. 1392E4 00 


I 
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TABLE  3 . c  RDX  MELT  W/OUT  BUBBLES 


1  H 

o 

300  K,  p  =  10 

MPa,  m  = 

30.2  kg/m2- 

s 

4+CDMDENSEE 

PHFlSE  PROFILES 

4-t 

J 

.  TEMP-OO 

I:  MLS  S'.  TEMP 

YNTR 

:  Y?  M ICR CMS 

KPT  -'TV  j  J.--' M£ 

1 

C.4785E+ 03 

0 .  0  0 

. 1  CODE 4 01 

0.  0 

0. 1E00E+ 0 

c. 

0. 4S94E4  03 

0.4  0 COE- 01  0 

.  1 CCCE+ 01 

0.  8 078 E 4 CC 

0.  1E62E4C 

•~r 

C.  ECC2E4  03 

0. 8CCCE- 01  0 

.  1  0  0  0E+  01 

0 . 4  076E4  0  C 

0.  1  688E4  IJ 

A 

H 

0.  31 1 1E4 03 

0.  1 £  0  0E4  0  0  0 

.  1  0  0  ,0  E  4  01 

0. 5999 E 4 CC 

0. 168EE4 0 

c 

0.  5219E+C3 

0.1 6 00 £4 00  0 

.  1 0  COE 4 01 

0. 78S2E+ CC 

0.  174  7E4  Ci 

if. 

0. 5328E4  03 

0. £CCCE4 00  0 

.  1 0  0CE+ 01 

0. 964 2E+ CC 

0.  1 8  09 E 4 0 

t' 

0. E437E4  03 

0.  £4  CCE4 00  0 

.  i 0CCE4 01 

0. 1 137E4  01 

0.  1871E4 0 

£ 

C.EE4SE4C3 

0. £8CCE4 00  C 

. 1 CCCE+ Cl 

0.  13C4E4  01 

0.  1932E+  Ci 

9 

0. 3834E-*-  03 

0- 3£ 0  CE4  0  0  0 

.  1  0  0  0E4  01 

0. 1466E+ 01 

0.  1994  E 4  0i 

1  C 

C. 576£E+ 03 

0. 36 CCE4 00  0 

.  1 CCCE+ 01 

0. 1624E+ 01 

0.  20E6E4  C" 

-  i 

j.  ± 

0. E871E4  03 

0. 4  0  0  0E4  00  0 

-  9S 9  0 f 

0 . 1 776E4 0 1 

0. 81 17E4  C: 

i  c 

0. 593 CE4 03 

0,44  00E4  00  0 

.  9 9 Si  E*t  0  C 

0 .  1 984 E+ 0 1 

0.2179E4  Ci 

0 .  6  C88E-+.C3 

0. 48CCE4 00  0 

. 9997E4 ct 

0 .  £ 068 E 4  0 1 

0.  884  0E4  Ci 

i  4 

0.  6  197E-+  03 

0.  38  00E4  00  0 

. 9994 £4 Ol 

0 . £2 08 E 4  0 1 

0 . £3 0 1 E4  0 

1  e 

X  _• 

0. 63CSE4 03 

C. S6CGE4 00  0 

. 9987E4 Cl 

0.234EE+C1 

0 .  £36 1 E4  Ci 

It. 

0. 6414E4 03 

0. 6CCCE4 00  0 

. 3976E4 00 

0.  £ 4  7 8  E  4  01 

0. 2419E4  Ci 

i  l' 

0. 65 33 E 4 03 

0-  64 0CE4  00  0 

. 99 3 EE 4 00 

0.26 08 E 4 01 

0.  247EE+  Ci 

i  3 

C. 663 IE* 03 

0. 68CCE4 00  0 

. 99 1 7E4  0  L 

C.272EE4 01 

0. 8E86E4  0: 

0. 674  CE-+  03 

0.7c  0  0  E  4  0  0  0 

. 98E1E4 CO 

0 . 886  Ct  +  0 1 

0.  85 68 E 4  O' 

£  C 

0.  ££4 8E-*  UJi 

0. 76CGE4 00  0 

. 97 26 E + CC 

0. 8983E4 01 

0.  2S26E4  Ci 

£  1 

0.  69E7E4  03 

0. 8CC0E4 00  0 

. 9E41E4 0 0 

0 . 3 1  CEE 4 0 1 

0 . 86  0 1 E4  0: 

in  c 

0. 7C66E4 03 

C.84  0CE+00  0 

. 92 1 EE4  0 C 

0. 3228 E 4 01 

0. 2E67E4  0i 

ill  w 

0 . 7 1 7 4  E4  03 

0. 88C0E4  00  0 

. 8678 E 4  0  C 

0. 33E4E+ 01 

Ci  .84  (  c'  E  4  ij  i 

dA 

0  ■  i  c.'c'-:'E+  Co* 

C.98C0E4CC  0 

.  781 CE4  CC 

0. 3487E4 01 

0. 2279E4  Ci 

“i  tr 

C.7391E4C3 

0. 96C0E4 00  0 

.  64  34  E 4  0  C 

0 . 3638 E 4 0 1 

0.  1938E4  Ci 

i££ 

0. 75CCE+ 03 

0. 1 CCCE4  01  0 

.431 3E4 CC 

0 . 3889E4  0 1 

0.  1379E4  Ci 

i 

HTDT,  J.--M3-S 
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NARROW  NEAP-FIELD  REGION: 
HIGH  GRADIENTS 

CONVECT I VE -D I FFUS I VE- REACT  1 VE 


Fig.  1.  Schematic  of  the  nitramine 
gaseous  flame  field,  showing  the  near¬ 
field  and  far-field  concepts. 


N1TRAMNE 
MASS  FRACTION 


Fig.  2.  Schematic  of  the  condensed 
phase  with  decomposition.  Gaseous 
bubbles  (decomposition  products) 
become  prominent  as  the  nitramine 
mass  fraction  drops. 
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APPENDIX  A 


CONDENSED-PHASE  DECOMPOSITION 
FOR  A  DEFLAGRATING  SOLID  PROPELLANT 
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Appendix  A 

Condensed-Phase  Decomposition 
for  a  Deflagrating  Solid  Propellant 


Two  distinct  cases  involving  condensed-phase  decomposition  of 
deflagrating  monopropellants  are  treated  herein  in  a  highly 
approximate  manner.  The  objective  is  to  obtain  simple  analytical 
expressions  for  the  extent  of  decomposition  (remaining  mass  fraction 
of  reactant)  and  the  thermal  gradient  at  the  surface.  In  the  first 
case,  assumed  to  represent  nitr amines,  the  extent  of  subsurface 
reaction  is  expected  to  be  small.  In  the  second  case,  assumed 
applicable  to  AP,  the  extent  of  subsurface  decomposition  expected 
is  relatively  large.  In  both  instances,  the  associated  activation 
energies  are  assumed  relatively  high,  and  the  process  exothermic 
overall. 

1 •  Ni tramine-like  Subsurface  Decomposition 

A  single,  first-order-overall  reaction  is  considered,  which  is 
consistent  with  the  observations  of  Robertson  [373  on  condensed- 
phase  decomposition  of  nitramines.  The  physical  model  and  coordinate 
system  are  depicted  in  Fig.  A.l. 

A  liquid  phase  near  the  surface  is  assumed,  where  all  subsurface 
decomposition  occurs.  Two  characteristic  features  are  assumed: 

(1)  only  a  small  extent  of  decomposition  -  up  to  15%  reactant 
depletion  within  the  pressure  range  1-10  MPa,  and  (2)  The  melting 
temperature,  T^,  is  sufficiently  smaller  than  typical  surface 
temperatures.  As  a  consequence  of  (2)  and  the  high  activation  energy, 
E  ,  most  of  the  reaction  would  occur  in  a  thin  region  within  the 
liquid  phase,  under  temperatures  close  to  T  .  Note  also  that  all 
products  of  decomposition  are  gaseous,  and  are  assumed  dissolved  in 
the  liquid. 

Let  the  dimensionless  thermal  enthalpy  and  the  length  coordinate 
be  defined, 


■C  sCT-T.VfTs-Tc^  Z=y/(v/mc=)  (A.l) 


The  associated  conservation  equations  for  the  region  of  the 
liquid  phase,  zm  <  z  <  0: 


where: 
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(A. 5) 


Kc=  /CcCTs-Tc) 

(A. 5) 

0  <  £STS/£C<<  X 

(A. 6) 

-KO/e  “£0 -fVfr+a,), 

1  (A. 7) 

and  Y(z)  denotes  the  reactant  mass  fraction. 


The  associated  boundary  data  are: 

Y(z»0  =  1 

(A. 8) 

,  rfo)-± 

(A. 9) 

■C-Cem  )  -T  'OZm  )  =  -<^m 

(A. 10) 

"hete  WW' 

The  last  condition  in  Eq.  (A. 10)  arises  due  to  energy  conservation 
at  the  liquid/solid  interface,  where  Q*,  the  heat  of  melting,  is 
depleted.  *n 

The  system  can  now  be  conveniently  transformed,  using  as  the 
independent  variable,  for  TT  <.  1: 


m 

Cr-<£)d<t>/ci-c  = 

(A. 11) 

Cc-4>)  dv/dT  =  hq  V 

(A. 12) 

where  cfiC 'C)=‘C— TT''  and  Y  =  Y  (tr) .  The  associated  boundary 

conditions  are 

VCcw)  =  1  (a. i3> 

As  stated  earlier,  most  of  the  chemical  reaction  is  expected 
to  occur  in  a  thin,  high-temperature  region  near  the  surface,  where 
ZT  ~  1.  Thus,,  the  quantity  0  <.£  <5sC  1  is  suggested  as  a  small 
perturbation  parameter,  and  the 't,— domain  is  conceived  to  have  a  narrow 
boundary-layer  (or  "inner")  region,  where  chemical  reaction  is 
important;  over  the  rest  of  the  liquid  phase  (or  "outer"  region) , 
the  reaction  is  negligible. 
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Sufficient  resolution  for  the  inner  region  is  obtained  by  use 
of  the  stretched  coordinate, 


17  =  (1  -r)/t 


Thus,  transforming  to  the  inner 
equations  are 


variables  <4?,  S'  and 


the  conservation 


_  -FCr})  — 

[l-<t>-sri)d(}>/dr\  =  ~£K*e  V  »-15> 
(i-4>-£ri)a|7«<r)  =  +  &K.€ K7|)  '»•“> 


where 


Feio- 


■7/a.i 

i-foya, 


A/ 


The  following  inner  expansions  are  suggested,  in  simple  powers 
of  £  : 


<W7J)'V  tb/ip  +■  &4^Cr))i-  (a. Hi 

f  ('»£>  =  iK-St-Jj)  ~  ^e0))+£|'(-r))+£2^C7l')  +  -  (A'181 

The  exponential  term  is  likewise  expanded  : 


(A. 19) 


Substitution  of  these  series  expansions  into  the  inner  conservation 
equations  (A. 15)  and  (A. 16)  results  in  an  ordered  hierarchy  of 
differential  systems  after  collecting  terms  in  equal  powers  of  £: 

Zeroth  order: 


=  0^  <%(<))  =  Co 

0  =0  ,  ^(0)  =0 


(A. 20a) 

(A. 20b) 


with  the  solutions 


n) = ^ = ^ AS<ri  %(*0 = 0 . 


ID 


(A. 21) 


Recall  that  according  to  assumption  (1)  herein ,  only  a  small 
extent  of  decomposition  is  expected.  Hence, the  boundary  condition 
for  "g~  is  homogeneous,  leading  to  a  trivial  zeroth  order  solution 
(implying  that  T  ^  1). 

First  order: 


— IQV^^W-C,  (A. 22a) 
=  ^  ,  <-CT)  =  |  (A. 22b) 


Note  that  now  both  boundary  conditions  are  inhomogeneous.  For  the 
reactant,  log  Y(0)  ^0(£)  is  expected;  accordingly,  a  defect  in  the 
incoming  heat  flux  from  the  gas  (relative  to  the  nonreacting, 
unperturbed  case)  must  also  be  anticipated,  owing  to  the  slight  extent 
of  exothermicity  in  the  condensed  phase. 

The  first  order  inner  solutions  are, 


t<ga.  -nta, 

t-  31  e 


K?a  i  1 


(A. 23a) 


+ 


£b,  -h 


Solutions  of  the  outer  region  where 
of  the  form 


Kcflt,  "]  (A. 23b) 

*n/  are  generally 


<=Kr)  ^4^Cr)+-  C  [M  Cr)  -e  e  <jg 


(A. 24) 


and  likewise  for  Y  (  ) .  These  are  similar  to  the  outer-field 

solutions  suggested  by  Bush  and  Fendell  £293;  detailed  calculations 
and  rigorous  matching  will  not  be  carried  out  herein.  However, 
important  insights  can  be  obtained  by  inspection.  Evidently,  the 
term 


^-fCO/e 


e 


for  typical  outer-field  values  of  *C  ,  cf .  Eq.  (A. 7) .  Thus,  the 
zeroth  order  solutions  are 


^tr)  =  ,  V0(.f)  =  1 
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determined  from  Eqs.  (A. 2) ,  (A. 3)  and  (A. 13),  using  the  outer 
expansions.  The  first  order  outer  solutions  would  obtain  exponential 
terms  analogous  to  those  in  Eqs.  (A.  23),  but  their  boundary  data  are 
homogeneous.  Consequently,  the  square-bracketed  constant  terms 
in  Eqs.  (A. 23)  for  the  inner  region  have  no  counterparts  in  the 
outer  solutions.  This  occurs  since  the  small  perturbation  (due  to 
reaction  near  the  surface)  can  not  penetrate  to  the  (outer) 
liquid/solid  boundary. 

Therefore,  the  present  approximation  can  be  concluded  by  imposing 


c,  =  Kc*d,/0-<fc’c')  -t-Ote)  <*-26> 

y,  =-  Kcffl./O-^O  +  OCe)  ,A-27> 

while  *  <pQ  =  -q^,  all  of  which  arise  due  to  inner-outer  matching 
requirements. 

The  important  consequences  of  Eqs.  (A. 26),  (A. 27)  are  that  the 
values  of  the  mass  fraction  of  remaining  reactant  G  and  the 
dimensionless  thermal  gradient  ~c'  (0)  can  be  defined  at  the  surface: 


<gCo)  =  %Lo)+  £f,C0)  =  £b<  =  -  £  +  OCt‘) 


<S>Crs,  m,T0)=  exp[-  £  ' 

Similarly, 


] 


(A. 28) 


r'03)=  l-4>Co)=  i-(4^e<t><(o))-i+^-e.^  +OO0 

—x  2 

Now,  since  e  *  1-x  +  0(x  )  for  x  «  1,  the  last  equality  may  be 

rewritten  in  terms  of  G: 


t'Co)  =  0+ck)- (i-s)siVQCn-Tc) 


The  explicit  expression  obtained  for  G  enables  now  calculations  of 
the  burning  rate  (using  the  eigenvalue  equation) ,  as  well  as  n(p) 
and  <5  (p) ,  where  partial  derivatives  of  ln(G)  with  respect  to  m, 
T and  are  required. 
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(NITRAMINE)  MASS  FRACTION 


SOLID 


'LIQUID  1' 


0 

Fig.  A.  1  Schematic  diagrams  of  condensed-phase  processes  described  herein 
for  nitramine-like  and  AP-like  monopropellants.  In  both  cases  a  thin  re¬ 
action  region  near  the  surface  prevails  due  to  the  high  activation  energie 


(AMMONIUM  PERCHLORATE)  MASS  FRACTION 


THIN 

LIQUID 

LAVER,  TC~TS 
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2.  Ammonium  Perchlorate-like  Subsurface  Decomposition 


The  formulation  herein  is  quite  similar  to  that  of  the  foregoing 
section,  and  the  same  general  assumptions  are  employed  regarding 
uriiformity  of  thermophysical  properties  and  the  high  activation  energy. 

The  following  are  exceptions:  (1)  The  subsurface  reaction  is  assumed 
of  order  zero  with  respect  to  its  dependence  upon  the  reactant  concen- 
ration;  this  implies  that  the  energy  equation  can  be  solved  independently. 
(2)  Melting  and  decomposition  are  intermingled,  so  that  the  "critical" 
temperature  T  (Tg)  denotes  herein  both  melting  and  the  onset  of  reaction. 
Consequently  (due  to  the  high  activation  energy),  one  may  expect 


C^-Tsy(T5-TO  «  i 

(3)  An  appreciable  extent  of  subsurface  decomposition  is  expected,  in 
contrast  to  the  previous  analysis  for  a • nitramine-like  substance. 

A  similar  analysis  was  carried  out  by  Cohen  [A.l]  and  previously 
by  Ramouhal  and  Cohen  [A.2]. 

The  energy  equation  for  X*c  <  ~C  <  1_  is 

Cr-4?)  Atp/dr:  =  (A.30) 

where  -Q  =  -Tc')/('  -q  -To  )  -  OC  O 

The  stretching  transform  of  Eq.  (A. 14)  can  be  applied  again  to  the 
entire  region;  after  neglecting  terms  of  0(£) ,  one  obtains  for 

o  <  •’is  Cl -•£)/£  <^lc 

the  energy  and  species  equations: 


Ci-cb)a4>/di^  =  -«*e  (a.3i) 

(I -<*0  V(71j=  1  (A. 32) 

which  obtain  the  solution 

<Kn)=  i-  [ 1  -2  Ka, 

(A.  35) 
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Note  that  only  the  negative  sign  of  the  square  root  term  has  been 
retained  [in  the  quadratic  equation  arising  from  solution  of  Eg. (A. 31)] 
since  at  the  surface,  <£(0)  is  expected  to  be  maximal. 

The  reference  species  (reactant)  solution  is  similar 

"rTTp^  C'4?C^)CcCTi-rT;V62t  (a. 34 a) 

where  C  is  a  linear  coupling  constant,  evaluated  at  the  "cold"  boundary: 

C  =  1  fc-TcYQ*  =  1  -Q*mfo£  (A.54>) 

Now  the  values  of  the  thermal  gradient  and  the  reactant  mass 
fraction  at  the  surface  n=0  can  be  evaluated: 

f  =t'co)  =  1  -  <kco= [1  +^+^„)+2«c*a,  {e 1)] 1/2 

CA.55) 

&  =  (1  -  QZ/Q?)  -4>io)<u  (Ti-T.y&f  = 

-  [(2*-Q*-ciCT5-T=)+-G:Ct^Tc')F  J/q*  (A.3&) 

where  G  is  expected  to  depart  appreciably  from  unity;  for  AP,  <5  == 0,3 
is  inferred.  Furthermore,  if  G=const  over  a  broad  enough  range  of 
T  (at  fixed  T  ) ,  then  in  this  range 

°  (36./3T5 \  =  0 

and,  by  Eq.  (A. 36),  this  condition  leads  to: 

$F/d ai  =  (i-fVCts-t0) 

F  ~  1  -  (Tf>-ToV0rs--T«O 

where  is  a  reference  surface  temperature,  and  the  differential  operator 
is  defined 

St/cLTg  S  ^/a rs  +  Cdm/dTs') 

All  of  these  partial  derivatives  can  be  calculated  from  Eq.  (A, 35), 

since  the  (T  ,  m) -dependences  of  each  term  are  known  explicitly,  cf.  Eqs . 

(A. 4)  -  (A.7f  and  (A. 10)*  Note  that  although  G  (T  ,  m)=const  for  fixed 

T  ,  one  may  still  obtain  for  variable  T  S 

o  o 
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The  foregoing  results  have  not  been  incorporated  in  the  calculation 
of  AP  temperature  sensitivity  in  the  present  study;  i.e.,  G=const  was 
assumed  under  all  conditions.  The  major  reason  is  the  uncertainty 
involving  T  ,  the  "critical"  temperature  (which  might  well  be  T  -dependent 
to  provide  5or  G=const) ;  without  a  clear  physical  definition  ofST  , 
calculations  would  tend  to  be  highly  speculative. 

It  should  be  emphasized  that  G (T  ,  m;  T  . . . ) =const=0 . 3  is  in  no  way 
inferred  in  this  study;  it  is  imposed^  according  to  indications  made  in 
the  analysis  of  Guirao  and  Williams  [32 J.  The  detailed  derivations  of  F 
and  G  in  this  section  would  be  more  useful  for  cases  where  G  is  variable, 
perhaps  in  simulation  of  the  so-called  foam  zone  in  double  base  propellant 
deflagration. 
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Appendix  B 

Sample  Calculations  and  Datum  Points  - 

Comparison  with  BDP  Monopropellant  Model 


1.  BDP  Calculations 

The  following  formulation  was  used: 

\  (b.d 

Zg  =  (B.2) 

yniyi')  =  <b.3) 

k;  =  (A/qOK,,  (b.4) 

f  =  (wV^eKf  (B.s) 

For  AP,  the  value  of  K°  was  obtained  by  fitting  a  single 
experimental  point  (at  p  =  4.5  MPa,  m  =  11.01  kg/m-s) ,  as  follows: 

T  =  884.7  k  from  Eq.  (B.3)  and  z  =  1.332  from  Eqs.  (B.l)  and  (B.2); 
tnus,  from  Eq.  (B.5), 

K*  (AfO  =  9.619  (ywvs)1/  ■Rx1'6’ 

Now,  according  to  Eq.  (B.4),  using  (  ^/C  )  =  5  x  10~5  kg/m-s,  on 
obtains  ^ 

1.924  g/cc-s-okJ'8 


G5 


which  is  the  kinetics  constant  for  gas  phase  reactjog,  only  slightly 
different  from  that  of  BDP,  namely  1.12  g/cc-s-atm  *  .  ^ 

For  HMX,  a  similar  procedure  (at  p  ■  3  MPa,  m  =  11.33  kg/m  -s) 
leads  to  T  =  1133  K,  Z  =  3.288,  and  hence 

S  B 

ICJ5 CHM*)  =  4.3*9 xIO-12  (l^/w^s)2/^ 


Eq. 


In  this  case, 
(B.  4) 


using  again 


( Vcp> 


5  X  10 


kg/m-s, 


one  finds  from 


(HM*)  *  <28678  ^/cc-s-ati^? 


2 

which  is  somewhat  larger  than  the  value  of  0.246  g/cc-s-atm  specified 
by  BDP. 


2 .  Present  Model  (MBR)  Calculations 


For  AP,  from  the  eigenvalue  equation, 


6  =  p'n'/ m/'  (B-6) 

=  (VCp/Ai/CE^Te/w)"1  <B-7> 

With  the  reference  temperature  ta^en  =  1000  K,  and  the  reference-point 
data  (p  =  4.5  MPa,  m  =  11.01  kg/m-s) ,  the  value  of  K  can  be  calculated 
from  Eq.  (B.6)  with  G  =  0.3  and  T  =  1400  K  (for  T  =873  K) : 

L  S 

=  I.387kW*"4 

Now,  from  Eq.  (B.7),  using  >/Cp  =  5  x  10  ^  kg/m-s  and  W  =  30  g/mol, 

A°C*P)  =  (.737 10  C^/w 


oC 


which  is  within  an  order  of  magnitude  of  the  original  Longwell-Weiss 
prefactor  when  converted  to  the  proper  units.  The  foregoing  value  of 
K^was  used  in  the  calculations  of  mjetc. 


For  HMX  the  strategy  is  slightly  different  and  aimed  to  obtain 
first  a  value  for  n^,  the  overall  gaseous  reaction  order  in  the  near 
field.  The  burning  rate  pressure  sensitivity  is  written 


^Gp)  *  OS  [  1  +  V(/l+-  B/p)]  ».») 


readily  obtained  from  the  HMX  burning  rate  correlation  31,38 

vn  =cf>p l4-e>/p3°-S 

with  B  =  3.84  fcJPa.  This  yields  at  the  reference  point  (p  *  3  MPa, 
m  =  11.33  kg/m  -s)  n  =  0.7193.  From  the  pyrolysis  formula,  T  =  683. 3K 
and  from  Egs.  (28),  G  (Tg,  m;  Tq  =  300  K)  =  0.9265  and  lg(G)  =  -0.0764. 
The  HMX  pressure  sensitivity  formula  is,  explicitly: 


cp>k|X 


Ts  Tl-'ZTo 
*  T^-Tc 


using  pertinent  data  from  Table  2?  thus 


2)jt^>  +  zJ 


(B.  9) 


hmx")=  \ao 

This  result  is  .in  remarkable  agreement  with  the  computed  value  of 
n1~  1.3  in  Ref.  38,  where  the  upward  shift  from  unity  is  explained 
by  thermal  enhancement,  due  to  the  presence  of  secondary  reactions. 
Now,  since  T^  =  constant  in  the  present  approximation,  the  eigenvalue 
equation  obtains  the  constant  factor. 


k°(HMX)  =  <SAriVpr?1  ■=  \.c\7KK)  7  ( 

evaluated  gt  the  aforementioned  reference  point,  with  n  =  1.4;  this 
value  of  is  used  in  the  m  vs  p  calculations  herein.  By  definition, 

<  S  [0a;  /( ]e  1 VTf 


G7 


where  values  of  particular  terms  In  the  square  brackets  are  unimportant. 
However,  to  facilitate  comparison  with  the  kinetics  data  (first  order 
overall)  of  Shaw  and  Walker  [”39]  ,  the  mean  temperature  is  taken  at 
T  =  1000  K,  and  a  multiplicative  factor  of  ^  is  included  to  adjust 
tne  overall  order;  the  prefactor  is 


1 2.7Sx  I0,sr 


k!?  ,  n,  etc.  This  value  is  with 
1/s  specified  by  Shaw  and  Walker. 


This  value  is  within  an  order 


with  the  foregoing  values  o|g  K 
of  magnitude  from  2.5  x  10  1/ 


In  summary,  each  of  the  models  considered  contains  one  overall 
kinetics  parameter  (or  two,  in  the  case  of  the  present  model  for  HMX) , 
which  was  adjusted  at  a  single  experimental  reference  point  (p,  m) . 
This  adjustment  is  specific  to  the  monopropellant  simulated.  For  the 
present  model,  the  adjusted  parameters  were  shown  to  obtain  kinetics 
prefactors  in  fair  agreement  with  available  chemical  kinetics  data. 
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NOMENCLATURE 


A  =  pre-exponential  factor,  Arrhenius 

kinetics  constant;  1/s  for  first  order, 
m  /mol-s  for  second  order 

C  ,Cc  =  specific  heats  for  gas  (isobaric)  and 
condensed  phase,  respectively,  J/kg-K 

D  =  diffusion  coefficient,  gas  mixture, 

m  /s 

E  =  activation  energy,  Arrhenius  kinetics 

constant,  kcal/mol 

h  =  specific  thermal  enthalpy,  J/kg 

k^/k^  =  reaction  rate  constant  (Arrhenius  form 

intended)  for  firs|  order  1/sec,  and 
for  second  order  m  /mol-s 

M  =  total  number  of  overall  reactions  in 

nonlinear  gas  phase  model 

2 

m  =  mass  burning  rate,  kg/m  -s 

N  =  total  number  of  chemical  species  in 

gas  phase  model 

n  =  burning  rate  pressure  exponent  in  apn 

formula 

p  =  pressure,  MPa  (Mega  Pascal; 

1  MPa  «  10.013  atm) 

Q  =  heat  of  reaction,  J/mol  (positive  for 

exothermicity) 

r  =  linear  regression  rate  of  burning 

propellant,  m/s 

<3S  =  surface  heat  balance  parameter 

(representing  liquid  side) 

Ru  =  universal  gas  constant,  1.987  cal/mol-K. 

In  equation  of  state,  MKS  units  used. 

T  =  temperature,  K 
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f 

G,  g 
i 

j 


=  final  flame  conditions 


gas  phase  property 


=  primary  decomposition  (i  =  1)  and 
secondary  reactions  (i  *  2,...M) 

=  chemical  species  (for  j  «  1,3... N-l) 
and  thermal  enthalpy  for  j  =  N 


L  =  liquid  phase 

m  =  property  corresponding  to  melting 

max  =  maximal 

s  *  liquid-gas  interface 

sub  =  sublimation 


v  =  vaporization  property 

0  =  ambient  conditions,  e.g.,  initial  solid 

phase  temp.  T^. 


Superscripts 


(  ) 


dimensionless  property 


(  ) 


*  mean  property 

=  specific  (per  unit  mass);  in  the  melt 
layer  analysis  only:  intrinsic  phase 
property. 

=  on  the  positive  (or  progressive 
coordinate)  side  of  an  interface 

=  equilibrium 
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u  =  mean  velocity  of  gas  mixture, 

perpendicular  to  propellant  surface, 
m/s 

W  =  molecular  weight,  kg/mol 

V*2  =  primary  an<3  secondary  reaction  rates, 

respectively,  mol/m  -s 

x  =  gY/  *,  reduced  liquid  phase  density 

fraction  in  melt  layer 

Y  =  mass  fraction 


y  =  length  coordinate,  m 


z  =  dimensionless  melt  layer  coordinate, 

Eq.  (12) 

^  =  dimensionless  near  field  coordinate 

in  the  gas,  defined  by  Eq.  (22) 


n 


*  log(^+  C),  dimensionless  gas  phase 
transform  coordinate 

=  thermal  conductivity,  J/m-s-K 


A  =  flame  speed  eigenvalue,  based  on  maximal 

primary  reaction  rate 


,co  "  '  i 

-  «V  V 


=  net  stoichiometric  coefficient  of  j-th 
species  in  the  i-th  reaction;  single 
and  double  primes  denote  reactant  and 
product,  respectively 


? 


«  mixture  density, 


kg/m 


3 


7^  =  (T-T  )/(T  -T  )  dimensionless  melt  layer 

s  m  2 

enthalpy 


=  dependent  variable  in  melt  layer 
formulation,  Eq.  (16) ?  also,  in  gas 
phase,  Eq.  (26) 


*  V  /V,  porosity  in  the  bubbly  melt  layer 


Subscripts 


B  =  bubbly  melt  layer,  overall  properties 

c  =  condensed-phase  property 
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